NUSC  Technical  Document  8559 
8  January  1990 


AD-A222  089 


A  Preliminary  Analysis  of  Threshold  Signal 
Detection  in  Ambient  and  Signal-Dependent 
Noise  Environments 


Weapon  Systems  Technology  &  Assessment  Department 


DTIC 


ELECTE 
MAY  3 11990 


Naval  Underwater  Systems  Center 

Newport,  Rhode  Island/New  London,  Connecticut 


Approved  for  public  release;  distribution  is  unlimited. 


Preface 


This  document  was  prepared  for  the  Weapon  Systems  Technology  &  Assessment  Department 
by  Dr.  David  Middleton  of  Sonalysts,  Inc.,  under  Contract  N66604-87-D-0546-0031. 


Reviewed  and  Approved:  8  January  1990 


F.  L.  White 

Head,  Weapon  Systems  Technology  &  Assessment  Department 


U.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


la.  SECURITY  CLASSIFICATION  AUTHORITY 


lb.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
TD  8559 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Naval  Underwater  Systems 

Center 


6c  AOORESS  (Cry.  sun,  *nd  ZIP  Cod*). 
Newport  Laboratory 
Newport,  RI  02188 


6b.  OFFICE  SYMBOL  1  7a.  NAME  OF  MONITORING  ORGANIZATION 

(If  app kc*bf*)  S 


Code  8219 


7b.  AOORESS  <OTy,  Sfaw,  and  ZIP  Cod*) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


Be.  ADDRESS  (City.  Star*,  and  ZIP  Cod*) 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  tpplkobl*) 

Code  8219 


1  10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO 

1 1  TITLE  (Inclua*  Socunty  CUssifKtaon) 

A  Preliminary  Analysis  of  Threshold  Signal  Detection  in  Ambient  and  Signal -Dependent 


11.  PERSONAL  AUTHOR(S) 


13a.  TYPE  OF  REPORT 
Tech.  Doc. 


David  Middleton 


Ub.  time  COVEREO  Il4.  DATE  OF  REPORT  (Y**r.  Month.  D*y)  hS  PAGE  COUNT 

FROM  7/1988  TO 2/1989  |  19913  January  8  I 


16.  SUPPLEMENTARY  NOTATION 


COSATI  COOES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  (Continu*  on  reverse  if  ntmury  *nd  lOtnnfy  by  bloat  numb*r ) 
Threshold  Signal  Detection;  Signal-Dependent  Noise, 
Reverberation,  Scattering,  Target  Modeling,  Doppler  Effect 
Wave  Surface  Scatter,  Optimum  Performance 


19  ABSTRACT  (Continu*  on  reverie  if  n*c*sury  and  idontify  by  block  numb*r) 

— Weak  signal  detection  of  signals  scattered  by  ocean  wave  surfaces  and  reflecting 
bodies  is  examined.  Optimum  threshold  algorithms  a  as  extensions  of  earlier  work  [l], 
[2],  [5],  [10],,'*are  determined,  which  contain  in  addition  to  the  usual  ambient  back¬ 
ground  noise,  signal-dependent  noise  generated  in  the  scattering  process.  Threshold 
performance  probabilities  are  formulated,  in  terms  of  the  now-generalized  performance 
parameters  associated  with  coherent  and  incoherent  detection.  As  before,  a  canonical 
theory  is  presented,  which  is  invariant  of  specific  noise  statistics  and  signal 
waveforms.  Both  spatial  and  temporal  processing  are  included,  the  former  implicitly 
in  many  cases  by  means  of  projected  beam  patterns. 

The  physico-geometrical  factors  of  platform  motion,  beam  patterns,  distributed 
spattering elements  and  source  and  receiver  geometries  for  both  the  "monos  taticr"  and 
"bistatic"'  regimes  are  summarized.  Various  distributed  target  models  are  presented, 
including  quasi-phenomenological  ones  capable  of  analytic  evaluation,  involving  both' 


20  distribution  /availability  of  abstract 
□  UNCLASSIFIED/UNLIMITEO  B  SAME  AS  RPT.  □  OTIC  USERS 


22*  NAME  OF  RESPONSIBLE  INOIVIOUAL 
Mr.  James  Kelly 


00  FORM  1473, 84  MAR  83  APR  edition  may  0«  used  until  e*heu*ted.  SECURITY  CLA 

All  other  edition*  «r*  obsolete.  Unciassir le 


21.  ABSTRACT  SECURITY  CLASSIFICATION 
Unclassified 


22c.  OFFICE  SYMBOL 

Code  8219 


PTJTf™ 

EmiwQki 


secumrv  classification  o*  tmi s  P»ce 

Block  19,  continued: 

"continuous  and  discrete  scatterers.  Under  tne  present  condition  of  narrow-band 
signals  and  far-field  geometries,  the  inpu|  signal  waveforms  are  explicitly  separable 
from  these  phys^co-geometric  factors.  A  ^weak^target  is  assumed,  in~<veak 
reverberation ,"H  vis-3-vis  strong  ambient  noise. 

Sufficiently  detailed  analytic  structure  is  provided^  along  with  the  availability 
of  ISD.^ro  achieve  explicit,  quantitative  results,  and  to  guide  the  treatment  of 
special  problems.  A  short  list  of  possible  next  steps  is  included. 


Unclassified 


SCCUBITV  CLA»SIFlCA  TlCN  Or  r»l$  PACE 


A  Preliminary  Analysis  of  Threshold  Signal  Detection  in  Ambient  and 
Signal-Dependent  Noise  Environments 

Table  of  Contents 


1.  Introduction 

1.1  New  Features 

1 .2  The  Basic  Geometry 

1 . 3  Aims  and  Organization 

2.  Threshold  Signal  Detection  in  Ambient  and  Reverberatory  Noise 

2. 1  Threshold  Detection  Algorithms  [1] 

2.2  The  Decision  Situation:  Signal-Dependent  Noise 

2.3  Performance  Measures  and  Bias  Terms 

2.4  An  Important  Special  Case:  Performance  Beams 

2.5  Analytic  Next  Steps 

3 .  Scatter  Signal  Models  0(s) 

3.1  Preformed  Beams 

3 . 2  Arbitrary  Arrays — Individual  Sensor  Elements 

3.3  The  Mean  Scatter  Signal  <0^> 

3.4  <0W>2 

3.5  The  Scatter  Signal  Covariance  <0$0$> 

4.  Elementary  Models  of  Target  Scatter  Signals,  (KT) 

4. 1  Target  Signals — Model  I 

4.1- 1  Coherent  Moments  of  0W 

4.1- 2  The  Incoherent  Moments  <C 

4.2  Model  la:  An  Approximate  Model  I 

4.3  Model  II:  "Point-Scatter"  Model 

4.4  A  Resonance  Component 

5.  The  "Telecommunication  Model”:  T— »x,  T— »S— >x 


6.  Results,  Remarks,  and  Next  Steps 
6. 1  Principal  Results 
.  6.2  Next  Steps 
Appendix  A.  1 
References 


1 

2 

3 

4 

4 

5 
7 

10 

11 

12 

13 

14 
16 
18 
21 
22 
25 

27 

28 
30 
32 

36 

37 

38 
38 

38 

39 

40 
42 


List  of  Principal  Symbols:  (See  also  Appendix  A.l  and  [5]) 


A0.  £<> 

a$,  ao.  aon.  a0,  etc. 
SoT.  SopT.  SopR 

«H.  «I 

(Xx.QrQp 

r(*)  „  . 

D  J  -coh,  me 
B(0s)(t) 

B(ox)(0 
£,  £0 
fo 
Ft 
g(;} 

H0,Hi,H'2,  H'i, 

I*x.  »x'.  etc. 

J  ~  MN;  j  =  m,n 
k0.  k 
KoCOin 
/,  In ,  f 

lcd.lw 

*1.  *2.  /3. 

M- 

mi,  m2,  m3 
ni,  n2,  03 
nR,  nA 

n 

A 

n 

X,  v0p,  v0T.  vt,  etc. 
¥ 

Q)o 

R 

R(>>  Rop 

S(t),  S0(t)in,  St 

5 


peak  signal  amplitudes,  (2.1b),  (3.1a) 

input  signal-to-noise  (amplitude)  ratios;  (to  receiver);  to  medium,  a0, 

directional  vectors 
false  alarm  probability 

direct  and  scattered  wave  fields 

beam  patterns,  or  beam  pattern  projections 
detector  bias  terms 

waveform  factor,  due  to  wave  surface  scattering 
waveform  factor,  due  to  target  scattering 
epochs 

carrier  frequency 
characteristic  function 
threshold  detection  algorithms 

hypothesis  states,  in  detection 
unit  vectors 

space-time  (sample)  indexes 
(angular)  wave  numbers 
signal  correlation  function 
nonlinear  detector  transfer  function 
(nongaussian)  noise  statistics 
direction  of  cosines  (x) 
p/q:  ratio  of  a  priori  data  probabilities 
direction  of  cosines  (y) 
direction  of  cosines  (z) 
normalized  (received)  noise  processes 
field  renormalization  operator 
surface  normal 

wave  numbers  (associated  with  beam  patterns) 

(received)  nongauss  +  gauss  noise  intensities 
2pf0  =  angular  carrier  frequency 
receiving  array  operator 
(plane- wave)  reflection  coefficients 
input  signal  waveforms 
detection  performance  parameter 

shadowing  function 


ii 


X 

0 

6i  9m,ni  6n 

0(s),  0(t) 

VD,  Avd 
wi 

Wrt-P,  WRT-p,  wrt-s 

C 


=  t2-ti 
=  error  function 
=  normalized  received  signals 
=  normalized  received  scatter  "signals,"  target  signals 
=  doppler  velocities 
=  noise  probability  density 

=  wave  number  spectra,  target  and  wave  scattering  surface 
=  wave  surface  elevation 


in 


Fig.  1.1a 
Fig.  1.1b 

Fig.  3.1 


Fig.  4.1 

Fig.  4.1a 
Fig.  4.1b 


List  of  Figures 

"Bistatic"  Regime:  direct  and  multipaths  to  target;  surface  "forward”  scatter  (Vc= 0). 
"Monostatic"  Regime:  direct  and  multipaths  to  target;  surface  and  target  backscatter 
(Vc=0). 

Geometry  of  Transmitter  (at  Op),  Receiver  (at  Or),  in  Vr,  and  the  m^-sensor  in 
Vr,  in  relation  to  the  effective  scattering  surface  A,  with  analogous  geometry  for  the 
p^-sensor  in  Vp. 

Geometry  of  the  scattering  body,  yis-k-vis  the  source  (T)  and  receiver  (R),  cf.  Fig. 
3.1. 

Geometry  of  the  "flat,  tilted  target"  model,  Ap. 

Point-scatter  target  model. 


IV 


A  Preliminary  Analysis  of  Threshold 
Signal  Detection  in  Ambient  and 
Signal-Dependent  Noise  Environments 

by 

David  Middleton 


1.  Introduction 

This  document  extends  the  earlier  work  of  the  author  ([1],  and  refs,  therein)  for  the  space- 
time  extraction  of  weak  signals  in  general  nongaussian  noise  to  include  signal-dependent  noise, 
i.e.,  reverberation,  and  distributed  targets,  in  an  underwater  acoustic  environment.  Although  this 
is  a  venerable  problem,  analogous  to  radar  target  detection  in  ground  clutter,  for  example,  recent 
analysis,  incorporating  statistical-physical  models  of  general  nongaussian  ambient  noise 
mechanisms  [2]-[4]  and  detailed  scattering  treatments  which  contain  the  relevant  geometrical  and 
physical  features  [5],  [6],  can  now  be  applied  and  generalized  to  this  class  of  problem  where  a 
significant  component  of  the  interfering  noise  is  signal-dependent,  as  well  as  often  nongaussian. 

In  addition,  the  results  of  this  study  provide  a  framework  for  obtaining  potential 
improvement  of  these  active  systems  by  means  of  "signal  design,"  whereby  a  subclass  of  signal 
waveforms  (subject  to  suitable  constraints)  may  be  found  which  further  reduce  the  probabilities  of 
decision  error  in  system  performance  [7].  The  extension  of  these  results  in  turn  can  lead  to  a 
possible,  more  restrict'd,  class  of  "lew  probability  intercept"  (LPI)  signals,  dependent  upon  the 
physical  constraints  involved  [8]. 

As  before  [11,  [9]-(l  1],  except  in  specific  examples,  the  treatment  is  generally  canonical,  in 
the  sense  that  the  functional  form  of  the  results  is  invariant  of  specific  signal  and  noise  models  and 
statistics,  as  well  as  specific  physical  circumstances,  geometries,  etc.  Fortunately,  this  is  always 
possible  in  threshold-signal  theory  [10],  [1 1],  which  in  turn  permits  direct  application  to  a  very 
wide  spectrum  of  practical  situations.  The  particular  physical  environment,  of  course,  provides  the 
needed  quantitative  calibration  to  the  canonical  results  [1].  [2],  [9],  [12].  Finally,  although  the 
analytic  emphasis  here  is  on  optimum  algorithms,  practical  sub-  or  "near-optimum"  procedures  can 
be  readily  derived  from,  and  compared  with,  the  former  [8]— f  10].  It  is  these  which  are  usually  the 
economical  ones  to  implement 


1 


1.1  New  Features 

Most  of  the  elements  of  our  present  analysis  have  become  available  during  the  last  decade. 
What  is  the  principal  new  feature  here  is  their  combination  in  a  single  formalism  which  permits  the 
controlling  details  of  the  physical  environment  to  appear  in  a  realistic,  i.e.,  statistical-physical  way. 
This  will  become  explicitly  evident  in  succeeding  sections.  Also,  because  detailed  derivations  of 
many  of  the  components  of  our  treatment  are  available  in  the  recent  literature  ([1],  [3],  [9]-[  12]  and 
refs,  therein),  we  shall  not  provide  them  here,  other  then  by  reference.  What  are  these  needed 
elements?  They  are,  concisely: 

I.  Nongaussian  Noise  Backgrounds,  with  more  or  less  significant  gaussian  components; 

II.  Ambient  and  Signal-Dependent  Noise-,  ocean  ambient  noise  often  has  significant 
nongaussian  components,  e.g.,  biological  sources,  shipping,  and  structural  noise, 
produced,  for  example,  by  arctic  ice  [13]— £15],  in  the  com  s,e  of  its  local  motion. 

III.  Distributed  Targets :  these  produce  signal-dependent  returns,  of  course,  which  consist 
of  both  resolvable  and  unresolvable  "multipath"  effects,  dependent  partially  on  the 
(partially)  random  orientation  of  the  reflecting  elements.  In  these  acoustic  cases  we 
may  also  expect  mechanical  resonance  effects,  excited  by  the  incident  radiation,  which 
can  be  a  significant  feature  of  the  target  structure.  Here  we  shall  introduce  some 
comparatively  simple  quasi-phenomenological  target  models,  whose  numerical 
evaluations  require  empirical  data  [cf.  Sec.  4.4  ff.]. 

IV.  Doppler  Effects:  these  arise  because  of  the  relative  motion  between  source  and  target, 
which  includes  both  deterministic  and  random  components.  Because  the  speeds 
involved  are  small  compared  to  the  (group)  velocity  (c0)  of  sound  in  the  ocean,  e.g., 
IvdI  «  c0,  a  manageable  analytic  theory  is  possible  [17].  For  any  realistic  treatment 
doppler  effects  must  be  included  here. 

V.  Explicit  Geometries:  quantitatively,  source,  target,  and  boundary  geometries  play  a 
critical  role  in  practical  signal  detection  and  design  problems  [1],  [5],  [8],  particularly  in 
our  present  class  of  problem,  as  is  analytically  evident  here,  cf.  Secs.  3.5  ff. 

VI.  Physical  Models:  these  also  play  critical  roles  in  our  analyses  since  they  determine  the 
(generally)  statistical  description  of  the  various  ambient  and  scattering  noise  and  signal 
mechanisms  in  the  course  of  propagation  in  the  medium  ([1],  [2],  [5],  for  example). 
Moving  wind-wave  surfaces,  ocean  bottom  elevations,  distributed  targets,  etc.,  are 
essential  mechanisms  here. 

All  the  above  have  received  detailed  attention,  recently  and  earlier  [cf.  [1],  [5],  and  remarks, 
refs.].  Our  principal  task,  as  noted  earlier,  is  to  combine  the  relevant  elements  of  I-VI  above,  to 
obtain  the  "macro-algorithms,"  so  to  speak,  and  their  expected  performance  measures,  viz.  here, 
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probabilities  r r  correct  decisions,  in  such  a  way  as  to  incorporate  the  controlling  physical  factors 
which  the  acoustical  environment  here  imposes. 

12  The  Basic  Geometry 

Figure  1.1  shows  in  a  symbolic  way  the  governing  geometry  for  the  class  of  problems 
examined  in  this  study.  A  general  "bi-static"  geometry  is  sketched  in  Fig.  1.1a,  while  Fig.  1.1b 
shows  a  typical  "monostatic"  configuration.  The  directional  lines  represent  effective  directions  of 
both  coherent  and  incoherent,  direct,  and  scattered  radiation  ([1],  Fig.  2).  For  the  present  we  shall 
ignore  volume  scatter  and  the  effects  of  velocity  gradients  (Vc  *  0)  in  the  medium.  This  restricts  us 
to  comparatively  short  ranges.  Also,  we  shall  assume  configurations  where  bottom  scatter,  if  any, 
is  quite  separable  from  target  and  wave  surface  returns,  and  can  therefore  be  disregarded  here.  All 
these  features  can,  of  course,  be  included  in  the  same  way  we  approach  the  problems  conforming 
to  the  geometries  of  Fig.  1.1:  the  results  are  simply  formally,  not  conceptually,  more  complex.* 


Fig.  1.1a:  "Bistatic"  Regime:  direct 
path  and  multipaths  to  targets;  surface 
"forward"  scatter  (Vc  =  0). 


Figure  1.1b:  "Monostatic"  Regime: 
direct  path  and  multipaths  to  target; 
surface  and  target  backscatter  (Vc  = 
0). 

H®T 


* 


This  is  not  entirely  true  for  the  longer-range  cases  where  Vc  *  0  must  be  accented  for.  See  [18]. 
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13  Aims  and  Organization 

Our  aims  here  may  be  concisely  stated:  they  are 

A.  To  obtain  threshold  detection  algorithms  for  scattered  signal  returns  in  signal-dependent 
(reverberatory)  noise,  as  well  as  ambient  noise  backgrounds,  when  the  latter  are  strong 
vis-a-vis  the  former, 

B .  To  predict  the  performance  of  these  threshold  algorithms,  in  terms  of  probabilities  of 
correct  decisions;  and 

C.  From  the  results  of  A  and  B,  to  indicate  an  approach  to  signal  design,  with  appropriate 
physical  constraints. 

As  stated  earlier,  the  analysis  is  directed  to  optimum  systems  (in  the  Bayesian  sense.  Chapter  1 8  of 
[7]),  whose  structures  are  to  be  approximated  by  practical,  or  "near-optimum"  algorithms,  which 
in  turn  are  usually  the  candidates  for  implementation. 

The  Report  is  organized  as  follows: 

i.  Section  2  provides  a  concise  summary  of  the  threshold  detection  formalism  extended 
now  to  include  signal-dependent  noise; 

ii.  Section  3  assembles  various  signal-dependent  (reverberation)  noise  results,  appropraite 
to  the  scenarios  of  Fig.  1.1,  while 

iii.  Section  4  introduces  some  quasi-phenomenological  target  models,  whose  distributed 
nature  is  specifically  noted,  along  with  a  resonant  component. 

iv.  Secton  5  contains  a  short  summai7  of  the  "telecommunication"  or  direct  path  cases, 
where  the  target  is  capable  of  detecting  the  emitted  signal. 

v.  In  Section  6,  the  role  of  performance  measures  is  discussed,  including  the  relevent 
physical  parameters  of  the  detection  scenario. 

vi.  In  addition  in  Section  6,  the  Report  concludes  with  a  short  critique  of  the  approach  and 
results,  along  with  suggested  next  steps. 

For  technical  details  and  derivations  the  reader  is  referred  in  most  instances  to  the  indicated 
references.  Reference  [5],  in  particular,  will  be  needed  here  also,  for  many  definitions  and 
descriptions,  in  order  to  avoid  an  excessively  lengthy  treatment.  Reference  [5]  is  recommended  as 
a  companion  document  for  the  detailed  applications  of  this  study. 

2 .  Threshold  Signal  Detection  in  Ambient  and  Reverberatory  Noise: 

Canonical  Structures 

Here  we  summarize  some  of  the  principal  results  needed  in  out  current  detection  situation, 
which  include  space-time  operation.  We  begin  with 
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2.1  Threshold  Detection  Algorithms  [1] 

There  are  three  principal  forms  of  optimum  (and  suboptimum)  threshold  detection 
algorithms,  depending  on  the  receiver’s  a  priori  state  of  information.  These  are  (1],  Sec.  VII: 


/.  Coherent  Detection 

In  this  case,  signal  epoch  (£<,)  is  known  precisely  at  the  receiver,  so  that  appropriate  a  oss- 
correladon  with  the  received  signal  is  possible.  The  corresponding  (optimum)  algorithm  is 

gj(x)  =  Bj.coh  -  (2.1) 

m,  n 


where 


/m.n=3^ogwl(x)]Xj 


m.n’ 


*m,n  —  x0n  hm). 


(2.1a) 


andxm  =  [xmi,  ...,xmn,  ...,xmN]  is  the  (normalized)  data  sample  at  the  111th  sensor,  xj  =  Xj/V'F, 
with  x  =  [xm],  m  =  1, ...,  M  sensors;  here  J  =  MN  is  the  total  number  of  space-time  samples.  The 
quantity  0m,n  is  a  signal-to-noise  ratio,  also  normalized,  such  that 


9m, n  =  aM(tn)s(m)(tn)  =  <s^}>  =  1 ;  (all  n),  each  m.  (2.1  b) 

and  a^  =  A^/V2T,  where  =  peak  signal  value  =  A^  (tn)  and  \jr  =  sum  of  the  intensities 

of  the  ambient,  signal-independent  noise  accompanying  the  signal,  including  system  noise  in  the 
receiver.  Bj.*0h  is  a  bias  term ,  cf.  (2.15a,b)  ff.;  also,  (64),  (66a)  of  [1].  Here  <  >  denotes  a 

(statistical)  average  over  any  ramdom  parameters. 

In  our  present  formulation  we  postulate  independent  sampling.  [This  can  be  realistic  for 
time-sampling  and  not  so  closely  approximated  for  spatial  sampling.  Analytically,  very 

considerable  simplification  is  thereby  introduced,  see  [1 1],  m  Sec.  10,  especially.)  Thus,  the  pdf 
wj(x)  of  the  data  vector  x  reduces  to  JJwi(xj),  j  =  mn  =  1, ...,  J,  where  wj(x)  is  the  first-order 
pdf  of  the  space-time,  here  signal-iniependent  noise  process  x,  cf.  (2.1a)  above. 

As  we  see  presently,  0  can  include  a  signal-dependent  noise  term,  cf.  (2. 10)-(2. 12),  which 
includes,  as  expected,  the  reverberatoray  or  scatter  component  form  boundaries  and  volume 
inhomogeneities,  and  the  possible  target,  when  all  are  weak  vis-i-vis  the  signal  independent  noise 
contributions. 

The  important  class  of  suboptimum  detectors  associated  with  gj.*oh  ,  (2. 1 ),  here  are  those 
for  which  some  simple  approximation  h(x)  of  /(x),  (2. la),  is  employed,  for  example,  strong 
"clippers,"  when  highly  impulsive  (ambient)  noise  is  the  dominant  interference  ([10);  [19],  Sec.  9). 
Thus,  the  corresponding  form  of  threshold  detector  to  (2.1)  becomes 
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(2.2) 


SJ(*)coh  -  B  j ,coh  -  Lhm,n<0m,n>»  hm,n  -  h(x)!x=xm  „• 
m,n 

(The  same  bias  is  used  in  (2.1),  (2.2),  to  insure  that  as  h  -» /,  gj  -»  gj.) 


II.  Incoherent  Detection  Algorithms  [1 ] 

This  "incoherent"  case  is  distinguished  by  the  complete  lack  of  knowledge  of  signal  epoch 
(£o)  at  the  receiver,  so  that  <0>  =  <0(£o, ... )>£q,...  =  0.  The  appropriate  threshold  algorithm  here 
is  found  to  be  ((60),  [1]) 

*  *  1 

SjOOinc  =  ®J-inC  +TT  £  S[/  mV  m',n'  +  /  m.nSmm'Snn’l  <0m,n0m\n’>,  (2.3) 

mmton' 

where  now 

<0m.n0m',n'>  =  pgl^  .with  m=  agag>  /  Hg  , 

p  =  <sg>sg>> 

/•  =  2L| 

m,n  (jx  xm,n  * 

The  quantities  m  and  p  are  respectively  amplitude  and  (normalized)  waveform  correlation  functions 
associated  with  the  input  signal,  either  directly  or  functionally,  as  we  shall  see  below  in  Secs.  3, 4; 

is  an  input  signal- to-(ambient)  noise  ratio.  Generally  rh  contains  any  fading  effects,  while  p  is 
influenced  by  "doppler-smear,"  which  is  primarily  a  phase  (or  frequency)  modulation  mechanism, 
cf.  (63),  [1].  For  the  corresponding  class  of  simpler  (sub-optimum)  detectors  /(x)  is  again 
repleaced  by  h(x),  cf.  (2.2). 


III.  Composite  Detection  Algorithms  [9]-[ II] 

This  is  the  most  general  threshold  formulation,  which  takes  advantage  of  both  coherent  and 
incoherent  information  about  the  received  signal  and  noise.  In  fact,  by  derivation  it  is  simply  the 
sum  of  the  coherent  and  incoherent  algorithms  (minus  an  a  priori  bias  correction),  viz.: 

e(J)-comp  =  g(J)-coh  +  e7)-inc-1°gli’  (cf.  Eq.  (2.2c),  [9]),  (2.4) 

where  u  =  p/q  =  ratio  of  the  a  priori  probabilities  of  there  being  the  desired  signal  in  the  data  x.  In 
this  way  any  coherent  signal  structure  can  be  used,  in  addition  to  the  associated  incoherent 
structure,  further  to  enhance  the  performance  of  the  latter. 
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22  The  Decision  Situation:  Signal-Dependent  Noise 

Using  Section  in,  IV  of  [1]  let  us  describe  the  decision  situation  where  now  signal- 
dependent  noise  is  included.  Since  the  medium  involved  is  linear,  this  is  formally  done  in  a 
straightforward  fashion.  We  have,  for  the  two  hypothesis  stated  here,  in  these  binary  detection 
situations  at  the  receiver 


Hy.SQN: 

X  =  Nr  +  Na  +  S  +  Sgcat 

(2.5a) 

vs. 

H0:N: 

X  =  Nr  +  Na  +  Sgcat» 

(2.5b) 

where  Nr,  Na  are  respectively  independent  receiver  and  (received)  ambient  noise,  S  represents 
the  desired  signal,  while  Sgcat  is  received  signal-dependent  noise  not  associated  with  any  desired 
target.  For  the  regimes  of  Fig.  (1. 1)  we  see  that  for  the  desired  signal 


St  -  (St)direct  +  (Sx)multipath.  (2.6a) 

where  (S^multipath  has  different  components: 


Monostatic  (Fig.  (1.1b): 

[(St)multipath  =  Sx-»S-n-»R+  Sx-»x-»S  -*r]  monoi 

i  (2.6b) 

Bistatic  (Fig.  (101a): 

[(Sx)multipath  =  Sx-»S-»T— »R+  Sx-»x-»S-»R]bi>  ■* 

*  (2.6c) 

Similarly,  we  find  that  for  the  interface  (S)  and  volume  (V)  scatter  that 

Monostatic  (Fig.  (1.1b): 

(SScat =  Sx-»S-»R+  ^X-»R)monc>;  1 

(2.7a) 

> 

Bistatic  (Fig.  (1.1a): 

(Sgcat  =  Sx->S-»R+  Vx-»R)bi.  ■* 

(2.7b) 

Of  course,  the  contribution  of  each  term  in  (2.6),  (2.7)  is  different  because  of  the  different 
geometries  involved.  All  of  the  components  are  signal-dependent 

However,  not  all  target  and  scatter  terms  are  comparable.  Even  for  rather  general  geometries 
the  multiple  interactions  with  the  random  interface  (S)  considerably  reduces  the  effective  return  to 
the  receiver.  Accordingly,  in  the  first  order  we  can  disregard  (Sx)multipath-  Moreover,  since  the 
interaction  with  wave  surface  (S)  [or  bottom  (B),  not  considered  here]  is  normally  much  stronger 
than  from  the  volume  inhomogeneities  in  these  geometries  where  the  surface  is  at  ranges 
comparable  to  the  target  (t),  we  can  ignore  volume  effects,  so  that  (2.6b, c),  (2.7)  are  replaced  by 
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(St)muldpath  =  0;  Sscat  =  St-»S-*R 


(2.8a) 


and  (2.6b)  reduces  to 


(St)  -  (St)direct  • 

Consequently,  the  decision  situation  (2.5a,b)  also  reduces  to 


Hj:Sm : 

X  =  Nr  +  Na  +  (Sx)direa  +  ST->S->R 

vs. 

_ _ 

X  =  Nr  +  Na  +  St-»5?->r  . 

(2.8b) 


(2.9) 


The  ambient  noise  Na  and  receiver  noise  Nr  cannot  be  ignored:  they  arc  always  present  and  set  the 
ultimate  limit  on  threshold  performance,  whereas  ST  and  Sscat  vanish  as  Sin  — >  0  at  the  source  (T). 
Our  rather  detailed  passage  from  Eqs.  (2.5)  to  Eq.  (2.9)  is  justified  here  on  the  grounds  that  it  is 
necessary  to  separate  the  present,  quite  common  cases  embodied  in  (2.9)  from  the  other,  more 
complex  and  rarer  ones  which  at  some  later  point  it  will  be  necessary  to  consider. 

Next,  in  normalized  form,  consistent  with  the  results  of  Sec.  2. 1  above,  we  can  rewrite 
(2.7)  as 


Hy.  ( s2®n .)  -> 

g 

x  =  nR  +  nA  +  e^t  +  e$-»s-*R)  1 

vs. 

Hq:  (sjGn:)  ->  Hj: 

x  =  nR  +  nA  +  e(s^s^R),  J 

(2.10) 


with  GW'W  =  0(Sin),  e.g.,  0  is  linear  source  function  (or  functional)*  of  the  transmitted  signal  Sin. 
Here,  specifically, 


and 


x  =  X/Vv;  nR  =  NrA/ \j/ ;  iu  =  Na^/y;  \|/ s <n£>  +  <N|>  , 


s  =  Vv®- 


(2.10a) 

(2.10b) 


It  is  important  to  note  that,  in  essence,  the  scattered  return  is  treated  as  just  another  signal,  albeit  a 


*  This  follows  at  once  Grom  the  fact  that  the  medium  and  scatter  process  are  linear,  cf.  Eqs.  (14)  and  (18)  of  [1]. 
®  il-sh'it  ft  s  j  lt/M- 
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random  one,  generally.  The  formal  effect  of  this  is  to  allow  us  to  replace  6  by  +  (fsHn  all  our 
earlier  results  involving  signal-independent  noise,  where  now,  of  course,  the  appropriate  physical 
structure  of  0(T)  and  8(s)  must  be  incorporated.  Formally,  the  signal  (S(T))  or  noise  (Hi  vs.  Ho) 

is  replaced  by  an  artificial  binary  two-signal  situation,  Hj  vs.  Hj,  cf.  (2.10),  where  0i  =  0(s)and 
02  s  0(s)  +  e(x).  For  these  cases  of  weak  targets  and  weak  reverberation  in  strong  nongaussian 

noise,  the  threshold  algorithms  (2.1),  (2.3)  are  at  once  modified  according  to: 


I.  Coherent  Detection  : 


<  0m,n  > 


6m.n  >=  <  Cn  >+<  0<m!„  >•  H2-  vs  <  >  «  <  0m.n  >-<  9£„  >.  «'l :  (2. 1 1) 


II.  Incoherent  Detection: 

<  0min0m',n'  >  -*<  0jT)6j'T)  >  +  <  0^0^  >  +  <  0(jT)0^  >  +  <  0(j.x)0(jS)  >,  H'2 ,  VS. 
<05T)0jx)>-4<05s)0j?)>,H,1;  j  =  mn,  j'=  m'n' ,  (2.12) 

where  due  account  must  be  taken  of  the  coherent  and  incoherent  components  of  the  target  and 
scatter  returns.* 

The  detection  algorithms  (2.1),  (2.3)  are  likewise  modified  here,  as  a  consequence  of 
(2.11),  (2.12),  to 

8coii  =  BJ-coh-X,j<ej,>>’  (2.13a) 

C  *  Bl-inc  ♦  jrl %  ♦  I) «JT>  [<  W  >  *  <  9<i”e</)  >  +  >]•  (2.13b) 


In  order  to  obtain  the  "anatomy"  of  the  noise  and  signal  terms  in  (2.10)  we  must  refer  to  the 
propagation  models  (cf.  Secs.  Ill,  IV  of  [1]),  from  which  we  see  that  (in  sampled  form) 


nA  “  ^aA>  ® direct  ~  R(aH)directl  ®(T->S-»R)  ~  **aI 

where 

a^U-fur^H;  aH  =  M«(-GT);  fisM^Qj 


(2.14) 

(2.14a) 


A  A 

in  which  Moo  is  the  integral  Green's  function  operator,  Qs  is  the  surface  inhomogeneity  operator 

A 

and  Gj  is  the  (point)  source  function,  which  includes  the  transmitting  aperture  or  array.  Here  R  is 
the  receiving  aperture  or  array  operator  (cf.  Sec.  II,  also,  Sec.  IV,  B  of  [1  ]).  For  surface  scatter  we 
may  use  the  results  of  [1],  Section  V,  Eqs.  (34),  etc.,  which  give 


*  We  remark  that  it  is  possible  for  <0*jT)0j^>  to  be  nonvanishing  when  both  0(T)  and  0(s)  have  a  common  specular 
component,  even  though  signal  epoch  is  not  known  at  the  receiver;  see  Sec.  4.1-2  ff. 
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A  A  A 

in  satisfactory  approximation,  since  R,  Moo,  Q  are  all  specifiable  in  detail  [vide  [5]  and  various 
results  below],  nA,  etc.  are  explicitly  describable  in  terms  of  the  propagation  physics  and 
geometry.  Thus,  our  attention  needs  to  be  directed  to  the  appropriate  field  forms  needed  to 
determine  the  components  of  x,  (2.10).  This  we  shall  do  in  the  remaining  portions  of  this 
document. 


23  Performance  Measures  and  Bias  Terms 

Along  with  the  (optimum)  threshold  detection  algorithms  summarized  in  Sec.  2. 1  above,  we 
need  corresponding  performance  measures,  as  well  as  the  bias  terms  appearing  in  (2. 1 ),  (2.3),  (2.4). 
From  [10],  Eqs.  (4.3a),  (4.5a)  we  write  directly,  using  (2.1 1),  (2.12): 


B!-coh=IW--^X  X{<0j2)>2"<0j1)>2}’  j  =  MN;j  =  m,n  (2,! 5a) 


m=l  n=l 


and 

>2 - < efef  >2].  (2 15b) 


Furthermore,  we  have 


°*j2  =  -2Bj  =  varH()  g]. 


(]1  *  p/q;  q  =  1-p,  or  ^21  =  P2/P1, .  (2. 16) 


cf.  (68)  [1].  This  last  (without  p)  is  the  key  parameter  describing  binary  performance  which  here 
is  expressed  in  terms  of  the  usual  Neyman-Pearson  (N.P.)  criterion,  as  the  probability  of  correctly 
detection  the  desired  siganl,  viz.: 


1  +  0 


£o j 

V2 


(2.17a) 


with  the  false-alarm  probability 


log(K/[i) 

V*Ci  J 


p  =  p/q 


(2.17b) 


and  K  =  a  preset  threshold,  compatible  with  the  a  priori  desired  value  of  otp.  [For  suboptimum 
threshold  algorithms,  cf.  (2.2),  etc.,  is  replaced  by  a0j  formally  in  (2.17),  but  different  values 
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for  L@),  LW,  cf.  (2.18)  ff.,  are  required,  [10]: 


L(2)  .  <  /2>Hq  ;  L(4)  »  <  wi7w}')2>Ho  ,  <  >0  *  Jwi(xIHq)  ( )  dx. 


(2.18) 


In  the  present  extension  of  our  work  to  include  signal-dependent  noise,  we  may  now 
combine  (2.11),  (2.12)  with  (2.14),  (2.15)  in  (2.16),  to  get  directly  the  following  general  forms  of 
the  (threshold)  detection  performance  parameters  Cj : 


with  s  =  V\/6,  cf.  (2.10b)  (2.19) 


Depending  on  geometry  <0js)>  may  or  may  not  vanish,  cf.  (Sec.  4.1-2)  ff.  Note  that  L(2),  (2.18), 
is  a  statistic  of  the  signal-independent  noise  nR  +  nA  only,  as  in  L(4)  also.  Similarly,  for 
incoherent  reception  we  obtain 


-  i  £{(l«> -2L<2>*)sii.+ 2L<2>‘  }[<  sf  >e</>  >2  +2  <  ><  ej-’af  > 

+2{<  ej*>e<?>  > + <  oftyp  >}{<  >*<  tftyp  >} + 9(*>e(p  > + <  9<u9<;>  >)2  ]  (2.20) 


Again,  the  cross-correlations  may  not  vanish,  particularly  if  the  specular  scatter  is  at  all 

noticiable.  This  is,  of  course,  strongly  dependent  on  geometry,  as  we  shall  see  below  [cf.  Sec.  4.1- 
2].  From  (2. 1 1),  (2.20)  the  critical  r6le  of  the  target  return  is  evident. 


2.4  An  Important  Special  Case:  Preformed  Beams 

In  the  case  of  preformed  beams  where  the  detailed  array  structure  is  implicit  in  the  beam 
patterns  (as  in  the  treatment  given  in  [5]),  the  spatial  sampling  of  the  signal  and  noise  fields  by  the 
discrete  (m)  sensors  is  now  continuous,  cf.  Appendix  II,  [5].  As  a  result  the  summations  Em  no 
longer  appear  explicitly  in  our  threshold  algorithms  and  perfomance  measures,  cf.  (2.1)-(2.4), 
(2.14)-(2.20).  Thus  0j  -» 0n  here,  and  we  have  the  equivalent,  alternative  forms  of  our  earlier 
theory  [9],  [10]: 


Eg-  (2.1) 


sVWcoh  =  BV-coh  *  ^/n<0(n)>;  In  = 


/(x  n) 


n=  1 


(2.21) 
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Eq.  (2.3) 


gj(x)inc  “  BV-inc  +  TT  n’ +  ^n'Snn')  * 

Art  jj* 

[<«W>  +  <  en)0($>  +  <  e(n)0(n)>]*  (2-22) 

with  corresponding  modifications  of  (2.3a,b,c),  i.e.,  formally  dropping  "(m)’’  therein,  i.e.,  J  — > 

N,  j  -»  n,  etc. 

Similarly,  the  bias  terms  B  (2.15a, b)  and  performance  measures  (2.19),  (2.20)  reduce 
formally  to 

Eq.  (2.15a)  =  log  p  -  I  [<9W>2  +  2<  0«><eg)>]  (2.23) 

n=l 

Eq.  (2.15b)->:  B^-inc  *  log  p  -r  I  {  (lW  -  2L<2)2)  6nn-  +  2L(2)2  }  • 

°  nn' 

{ <0W>2 + 2<e5)ei|)><eW> 

+  2[<e(?e§)>  +  <  0<?eg)>]  [<eW0§?>  +  <  e^>eff>>] } 

(2.24) 

where,  by  (2.16),  o^2  and  are  related.  The  former,  by  (2.11),  (2.12),  j  -> n  here,  become  the 
results  below,  equivalent  to  (2.19),  (2.20),  viz.: 

Eq.  (2.19)  O$<oh  =  L<2>  £  {  <9<?>2  +  2<  0<J>><0$>  }  (2.25) 

n 

Eq.  (220)  a^.inc=j  £  { (L<4)  - 2L(2)2)  8nn-  +  2L(2)2  } 

{ <0W>2 + 2<o  W><eWe§)> 

+  2  [<e(s)0(n^><0«Q(s)>  ]  [<0(*)0§)>  +  <  e^>e5)>] } . 

(2.26) 


The  pertinent  results  of  Sections  3  and  4  ff.  are  applicable  here. 

2.5  Analytical  Next  Steps 

The  next  steps  to  be  taken  in  our  analysis,  which  are  required  to  give  explicit  structure  to 
both  the  bias  and  performance  parameters,  cf.  (2.16),  (2.19),  (2.20),  are  concisely  stated.  They 
are: 

I.  Obtain  0®  and  <0^>,  <0®>2,  for  coherent  (and  incoherent)  detection,  cf. 
(2.19):  ", Scatter  Signal  Models"  (Sec.  3) 
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II.  Obtain  (0®  and  <0®>,  <0(j)>2,  also  for  coherent  detection,  cf.  (2.14)  again: 
"Target  Models"  (Sec.  4) 

III.  Construct  <0<j¥ji>>,  <0<j¥j>>,  and  <0^j¥f)>,  cf.  (2.20),  needed  for 
incoherent  detection. 

Specifically  here  we  shall  make  the  following  assumptions: 

i .  "narrow  band"  input  signals; 

ii.  far  field  (i.e.,  Fraunhofer)  geometries; 

iii.  Kirchoff  theory,  and  associated  approximations ,  described  in  Sec.  II,  E  of  [1]. 

The  narrow-band  assumption  permits  analytical  solutions  for  the  "signal"  terms  0(s)-  0to  in  (2.19), 
(2.20).  (However,  at  a  later  stage  we  may  expect  to  have  to  consider  broad-band  cases,  where 
beam  patterns,  absorption,  etc.  are  frequency  dependent  and  where  final  results  must  be  obtained 
computationally.) 

Also,  for  the  LPI  problem,  we  shall  need  to  consider  the  "detection  recognition"  problem, 
where  the  target  (now  the  receiver)  does  or  does  not  detect  the  transmittal  signal.  This  formally  has 
the  effect  of  setting  0(s)  =  0  and  replacing  0(T)  in  (2.19),  (2.20),  etc.  by  0j,  the  received  signal  (at 
the  target)  from  the  original  source  (T).  As  we  shall  see,  this  is  the  simple  one-way  direct  path  or 
"telecommunication  model,"  which  does  not  involve  a  distributed  target  configuration  [Sec.  5  ff.],* 
although  it  does,  of  course,  require  an  appropriate  receiving  array,  or  beam,  cf.  Secs.  3, 4. 

Similarly,  for  the  important  case  of  preformed  beams  (Sec.  2.4)  we  parallel  I-III  above, 
now  for  0®  -»  0(j[),  0(p  — >  0W. 

3.  Single  Scatter  Models,  0(s> 

We  consdier  here  first  the  "classical"  cases  of  preformed  beams  where  the  details  of  the  array 
structure  are  implicit  in  the  resultant  beam  patterns  Qj  (v0  -  v0Tlfo)»  etc.,  and  correspondingly  in 
the  beam  pattern  projections  on  the  (mean)  scattering  surface  S0(x,y),  e.g.,  dp  -*  Qt  (r,ko),  etc., 
cf.  (17>— <19)  vs.  (20b),  (26)  of  [5],  where  QrtIs0  is  the  combined  beam  pattern  projection  of  Qt 
and  Or. 

In  this  formulation  0j  — *  0„,  cf.  Sec.  2.4  above,  but  the  spatial  decorrelation  of  the  scattered 
field  on  S0  can  be  extensive,  i.e..  At  ~  A  »  Xo2  ~  X<2>  where  Xc  is  an  effective  correlation 
distance  of  the  random  surface  £(r,t),  measured  on  S0,  and  A  is  the  effective  "illuminated-viewed" 
surface,  vide  Fig.  1  of  [5].  Thus,  the  spatial  processing  gain  inherent  in  spatially  independent 
samples  (m)  is  implicit  in  the  beam  pattern  projections  on  the  scattering  surface  S0(x,y),  when  A 
»  Xo2  ~  Xc2  so  that  many  "independent"  scatter  samples  are  available. 


*  If  surface  scatter  is  also  important,  because  of  source  and  target  geometry,  then  we  may  replace  0  by  0  +  0W,  in 
(2.11),  (2.12),  (2.19),  (2.20),  etc. 
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3.1  P  ref  on  .  »C  d  Beams 

From  (20a),  (26)  of  [5]  we  have  directly  for  the  received  scatter  signal  sample  in  the  scatter 
geometries  described  by  both  Figs.  (l.la,b)  and  Fig.  1,  [5],  when  multiple  scattering  (in  S0)  is 
ignored  (/  =  0): 


Eq.  (20a),  [5]: 


0<ns>  =  Re(S<*W^f) 


•c(s) 


Aq  ^  A o 


9(ns)  =  Re 


& 


=  V2^“n 

«»i+d(20)  . 

JSi„(s/2jd)F|0)(slC(r,tn),...)es(t»-£)^ 
— i+d(=Br,) 


tn  =  nAt. 


(3.1a) 


where  Sin(s/27ri)  is  the  amplitude  spectrum  of  the  normalized  input  signal  waveform,  to  the 
medium  and  where  A0  =  peak  value  of  that  signal  (which  is  different  from  A’0,  (2.16),  =  peak 
signal  input  to  receivers. 


Eq.  (26).  IS 7:  F§»  =  j ^„S  G+s  fl  RT(r.k0) 


with 


e-(s/co)[(r+0  •  2ao+cT0]  .  e-(s/c0)D^(vD/tn-e,r...)  dxdy> 

dS0  =  dydy;  ko  =  C0o/co  =  27cfo/c0;  r  =  ixx  +  iyy  on  S0;  T0  =  (Rot  +  RoR)/c0; 

Gd)  s  e‘2a®oc°T°/(4a)4R0y  +  Rq2r  =  Eq.  (A.  1*3) 


(3.1b) 

(3.1c) 


Here  D0  is  the  doppler  (displacement)  contribuution  from  the  platforms  (T,R),  surface  drift, 
deep  current,  etc.  In  the  present  model  D0  is  found  from  (24),  [17]  to  be  specifically 


D§  3  (2a0*vc)(tn-E)  +  (2a0*vd)(tn-£-RR/c0)  +  {-( i0rvoT)(tn-e-T0)  +  ( i0R*v0R)(Ve)  1(3.2) 
volume  drift  surface  drift  platform  motion 

=  "deep  current" 

where  vc,  vd,  v0t»  v0r  are  respectively  the  "deep"  or  volume  current  (if  any),  wind-induced 
surface  dnft,  and  platform  velocities.  Thus  Rr  =  Rqr  -  IqR  -r,  and  because  here  Irl  «RoR  and 
the  dopplers  are  small,  e.g.,  Ivl/c0  «  1,  we  can  let  Rr  ->  RoR  in  D0,  so  that D0  is  essentially 
independent  of  r  vis-k-vis  the  effects  of  the  moving  surface,  ~(r  +  z).  In  addition,  these  velocities 
may  have  random  as  well  as  deterministic  components.  The  velocities  themselves,  however,  are 
here  assumed  to  be  independent  of  both  space  and  itme  variations;  (these  can  be  introduced  later, 
when  needed). 

Also,  in  (3.1)  we  have  (cf.  Appendix  A.l  ff.) 
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20o  =  10t - i0R,  Eqs.  (A.l-3),  etc.  of  [5];  Gd>  =  Eq.  (29b),  [5];  (A.l-3)  here.  (3.2a) 


where  Fig.  3.1  below  shows  the  directional  unit  vectors  ioT.  »oR-  The  wave  surface  elevation 
(above  S0)  is  given  by  £(r,t)  as  before.  The  quantity  Sin  is  the  amplitude  spectrum  of  the 
normalized  input  signal  waveform  s^  ->  sn  =  s(tn),  cf.  (2.1b)  applied  at  each  element  (m)  of  the 
transmitting  array  implicitly  described  by  Ox.  The  first  exponential  term  in  (3.1b)  represents  the 
path  delay  and  angle  modulation  in  the  scattered  emissions,  caused  by  the  surface  motion,  while 
the  second  contains  the  effects  of  any  doppler  shift  and  spread,  embodied  in  vq. 

The  effects  of  both  transmitting  and  receiving  beam  patterns  are  contained  in  their  joint 
projection ,  QrtIs0»  on  So.  cf.  (3.1b),  which  may  include  many  correlation  distances  (~  X<)  of  the 
random  wind-wave  surface  about  S0,  cf.  the  introduction  to  Sec.  3  above.  Appendix  A.  1  provides 
definitions  of  the  various  other  elements  of  (3.1a,b). 


Fig.  3.1.  Geometry  of  transmitter  (at  Ot),  receiver  (at  Or),  in  Vr,  and  the  m^  sensor  in  Vr,  in 
relation  to  the  effective  scattering  surface  A,  with  analogous  geometry  for  the  p*  sensor  in  Vj. 

For  the  narrow  band  signals  used  here, 

Sin(s/27Ci)  =  Sin  0,  Is'l  «  s0,  (3.3) 

V27ti  2m ) 

which  allows  us  to  replace  s  by  s'  in  (3.1a,b),  with  es^n_e^  =  eso(trr€)+s(tn-€)  jn  (3.ia)  in  the  usual 
way,  viz.: 
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(3.4) 


9(?  Inn  =  Re  {  a^S0(tn')in  F<§)(iko,tn,...)eit0otn’  }  ,  t„* 
’  ff 


where  we  have  replaced  S  by  s  again  in  (3.1b)  also,  and  where  F(°)  represents  (3.1b)  now  with 
e~sT°  removed.  Here  S0(tn')  is  the  complex  envelope  of  the  (normalized)  driving  signal  s(tn), 
while  F(°)  contains  the  physical  structure  of  the  propagation  process.  We  now  note  the  critical 
feature  of  (3.4):  the  input  signal  and  the  physical  structure,  Fl°K  are  separable,  for  these  narrow- 
band, far -field  cases,  even  when  the  doppler  is  not  independent  of  position  [cf.  remarks  following 
(3.2)]. 

32  Arbitrary  Arrays — Individual  Sensor  Elements 

We  turn  next  to  the  important  cases  of  arbitrary  arrays,  considered  on  the  basis  of  hte 
individual  (m*)  sensor,  so  that  the  general  algorithmic  and  performance  results  of  Secs.  2. 1-2.3 
may  be  applied  directly,  when  the  0®,  etc.  are  explicitly  obtained. 

From  Eqs.  (30),  (31),  of  [1]  we  can  write  for  the  contribution  of  the  m*  receiving  sensor  to 
the  receiver's  input 

Or  -» IA^f0)l  .  e2,ti(XRm-yoR),rm>  (3.5) 

where  jv'or  (=  i0Rfo/c0)  is  the  (vector)  wave  number  representing  the  beam  steering,  as  earlier,  and 
where  the  directional  wave  number  of  an  elementary  scattering  source  at  r,  cf.  Fig.  (3.1)  is, 
namely, 

VRm  =  iRmfo/co>  (3.6a) 

fam  s  (RoR  +  rm  ~  **)  /  IRoR  +  rm  ~  *"1.  (3.6b) 

Using  the  far-field  condition  and  the  fact  that  the  array  size  (~  Lmax)  and  die  correlated  portions  of 
the  scattering  surface  A  are  small  vis-^-vis  Rot.  RoR  allows  us  to  write 

?Rm  s  t0R  -  (r-rm)  •  (and  1r  s  10r  -  r  •  similarly),  (3.7a) 

where 

SoR  =  I  -  «oR'oR  (3-7b) 


16 


is  the  symmetrical  diadic  (or  3  x  3  matrix) 


(*oR)x(^>R)y  (*oR)x(*oR)z 


a0R  = 


(•oR^'oR^y  i'(*oR^y  (*oR)y(*oR)z 
OoR)x(‘oR)z  (*oR  )y(*oR  )z  I-(‘oR)z 


whose  specific  elements  are  readily  obtained  from  (A.  1-2)  of  [5],  etc.,  for  general  and  monostadc 
configurations.  The  (complex)  weighting  A(m)  depends  only  on  the  signal's  center  frequency  (f0). 
Our  normalized  scatter  siganl  8®  now  becomes,  from  (3.1),  (3.4),  with  (3.5),  for  these  n.b. 

signals,  in  the  far-field  regime: 


now  with  (3.5). 


^(ik0|...)m=Wik”D”(VD,,--e'-).  aT(r.k0)A<f)(f0) 

J  V  nz  Jg+s 


.  L^o  (-i0R  -i'oR  +r’  <*oR  /RoR  )} •  rm  .  rm  -aoR  •  rra  /  RoR  1  c-ik0 (r+Q  •  2a0 dxdy; 


S-*  1C|>0 


(3.8b) 


The  exponential  term  ko{ }  is  the  brackets  []  of  (3.8b)  is  the  wave  number  associated  with  the 
steered  m*  sensor,  including  radiation  from  all  points  (r)  on  the  scattering  surface  (S0),  while  the 
second  exponential  contains  the  correction  for  the  displacement  of  the  m^1  sensor  from  the 
reference  Or,  cf.  Fig.  3.1.  The  explicit  relations  in  rm  show  how  the  sensor  locations  can  modify 
the  Total  Surface  Scatter  Function  (TSSF),  Ff°),  cf.  (3.1b),  (3.6b).  Beam  formation  here  is 
achieved  by  summing  over  (m),  e.g.,  Lm,  as  shown  in  (2.1),  2.3),  and  (2,19),  (2.20).*  Again, 


’  However,  it  is  clearly  not  true  the  0n  =  In,  n,  since  £m<0m  n>^  *<(£m6m  n)^>-  etc.,  vide  (2. 14),  (2.15), 
(2.29),  (2.30). 
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the  input  signal  and  the  physical  structure,  embodied  in  the  TSSF,  F(°)  here,  factor  explicitly,  cf. 
(3.4),  with  D0  given  by  (3.2)  in  these  cases. 

In  similar  fashion  we  can  "anatomize"  the  transmitting  beam  pattern  t.  and  write  specifically 

here 


dT  =  XIATP)(fo)|e"i^)(fo)'e2,ti(UTp‘U,oT)*rp  => 

P 


ax(r,k0). 


(3.9) 


this  last  for  the  projection  of  on  S0,  with  v'0j  =  i0Tfo/co-  the  steering  wave  number,  and 
following  (3.4),  (3.5) 


vTp  =  *Tpfo  / coi  *Tp  =  *oT-(r-rp-aoT/RoT:  a0T  =  I_loTioT/Cf.  (3.5c),  (3.9a) 


where  again  we  employ  the  far-field  and  narrowband  conditions.  Note  that  because  of  these 
conditions,  and  the  result  (19),  [5],  (X-\  necessarily  appears  as  a  beam  (£p)  and  thus  as  a  beam 
projection  on  So,  in  the  TSSFs  F^  and  F^m- 

The  importance  of  the  "anatomization"  of  r  and  t»  lies  in  (i)  application  of  specific  practical 
configurations  of  sensors,  and  (ii)  to  beam  optimization  by  direct  "matching"  of  the  sensor 
configuration  to  the  spatial  statistics  (covariance,  usually)  of  the  ambient  and  scattered  fields  ([8]; 
Sec.  10,  Sec.  10.4  of  [1];  [20],  [21].  This  incorporation  of  spatial  matched  filtering  with  its 
venerable  temporal  analogue  (see  Sec.  16.3  of  [7]),  is  currently  being  called  "matched  field" 
processing. 


3.3  The  Mean  Scatter  Signal  <  Q( S)> 

We  now  carry  out  the  statistical  averages  denoted  by  <  >,  viz.: 


<  > 


<>e;e'  =  [S’o.C.Cx.yVD,,  J  , 


(3.10) 


with  <  >e  =  J8(e-eo)de  and  where  we  use  the  assumptions  and  approximations  of  Sec.  3,  [5], 
namely  £x,  £y  are  independent,  as  are  all  the  doppler  components  in  (3.2),  as  well  as  ao,  etc. 
Then  (3.4)  becomes 


<9®>  =  R={yo(fon  >in  <F|°>(fon  } 


t  on  s  n^t  To.  (3.11) 
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where  now  specifically 


<  F|o)  >  =  ik0VG^>/^S(-2a0Z)elk°15«:(l'°")(F1)AVDF1(-k02a02)2  •  wRX(2k0ot01)s.  (3.1  la) 


The  various  components  of  (3.1b)  — >  (3-1  la)  are: 

=  “2cco2  =  b0  =  cos0OT  +  cos0OR,  Eq.  (4.2d)  and  (A.  1-3)  of  [5] 
CxCy:G+S=G  (3.12a) 

<e-ikoC*2Oo>£  =  Fi(-2koOoZ)c,  (3.12b) 

the  (first-order)  characteristic  function  (c.f)  of  wj(Q, 
the  1st  order  pdf  of  surface  elevation,  £,  which  is  here 
assumed  to  be  homogeneous,  i.e.,  independent  of  r; 

=  e-R(Gsgy2,  Rg  s  (2koCtozO;)2,  aC2  =  Q,  <>  =  0  (3. 12c) 

when  £  is  a  gauss ian  process,  which  is  a  good  model 
here  for  the  large-scale  ocean  wave  surface,  cf.  (39), 
(40)  of  [5].  Thus  surface  roughness,  as  is  well 
known,  acts  to  destroy  coherence  of  the  scatter  signal. 

e-ikoDoc(t0n)  =exp  [-ik0(2a0»vc)(t'on+T0),exp  [-iko{-(2(Xo*vd)(t’on+ToRoR/Co}] 

•  exp  [-ikof-aoT'v^Tlt'on+OoR’VoRjd’on+To)}]  (3.12d) 

which  are  the  mean  doppler  components  of  deep 
current  (vc),  surface  drift  (vd),  and  platform  motion 
v0T.  v0r: 


(FI)Avd  S  (Fi)Avc*  (Fi)Avd*  (Fi)AvqT  *  (Fi)Av0r,  (3.12e) 

wl.  re  each  (Fi)Av  is  the  c.f.  of  the  fluctuation,  Av,  in 
the  various  velocities  of  the  present  doppler  model. 
Specifically,  we  have  in  the  case  of  the  usual  gaussian 
pdfs  of  the  Av’s: 


'  (Fl)Avc  =  exp  {-ik2(t,on+To)2[(2a0xOcx)2  +  (2aoyaCy)2  +  (2ctoz<Jcz)2 ]/2)  ,  (3. 120 

j  (Fi)Avd  =exp  {-ik^t'on+To'rloR/Co^^OCox^dx)2  +  (2cCoy<j(jy)2  +  (2aozt7d z)2  ]/2}  ,  (3.l2g) 
1  (Fi)AvoT  =  exp  { -ik2(t'on)2[(2aoxaTx)2  +  (2otoyOTy)2  +  (2aozOTz)2  J/2 }  ,  (3. 12h) 

(Fi)Av0r  =  exp  {-ik2(t'on+T0)2[(2aoxORx)2  +  (2aoyORy)2  +  (2aoz<JRz)2 1/2}  ,  (3.12i) 

where  the  Cx.y^'s  here  aare  the  respective  variances  of 
the  doppler  components.  The  critical  feature  of  these 
doppler  terms  is  the  progressive  degradation  of  <0(s)> 

A 

as  the  duration  (T)  of  the  signal  (~S0(t'0n))  is 
increased,  a  degradation  produced  by  the  random 
doppler  fluctuations,  Av,  as  can  be  seen  directly  from 
(3. 12f-i)  above.  In  the  limit  T  °o,  (Fi)av  — »  0,  and 

as  expected  this  doppler  "smearing"  destroys 
coherence  completely,  e.g.,  <0(s)>  — >  0.  For  "short" 
signals,  of  course,  this  effect  can  be  small,  provided 
the  0's  are  themselves  sufficiently  small. 

wRT(2koa0j.)s  =  wRT(2ko0Cox,2kQaoy)s ■* J#RT(r)e2,koSo*rdr,  dr  =  dxdy.  (3.12j) 

So 

This  is  the  (amplitude)  wave  number  spectrum  of  the 
joint  beam  pattern  projections  of  the  reference  surface 
S0l  with  ctoi  the  surface  components  of  cto,  since  r 
has  no  z-component 

For  the  gaussian  beam  patterms  and  projections  used  in  our  current  model  (cf.  Appendix  II  and  Eq. 

(A.2-19)  et  seq.  of  [5],  where 

RT(r)  =  gTgRe'kor  •  BTRT/2  +  *0r-bTR>  (3. 12k) 

(3.12j)  becomes,  with  the  help  of  (A.2-13b),[5], 

wRT(2k0q0i)s  =  e-f(2aox-bTRx)2/A  +  (2aoy-bTRy)2/B}/2  (=  [!£  Eq.  (41)t  [5]) 

=  2gTgRAREFexp  {-^{(2(X0x-bTRx)2/A  +  (2aoy-bTRy)2/B]),  (3.121) 
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from  (A.3-7)  of  [5],  where  Btr  is  the  (dimensionless)  (2x2)  matrix 


Btr  = 


rA  °  i 

J  *  Eq.  (A.2-2)  of  [5]: 


A  =  At/RJt  +  Ar/R^,  B  =  a^-AT/R^x  +  ^Ar/R^r 


bxRx  =  C;  bjRy  =  D 


(3.12m) 


cf.  A.2-2b  of  [5]  for  details,  and  Eqs.  (A.2-19)-(A.2-22)  therein;  usually  bxRx.y  =  0.  From 
(3.121)  is  is  clear  that  unless  the  beam  orientations  are  essentially  in  the  specular  direction,  so  that 
2oq  -  bxR  =  0,  wrt  becomes  exponentially  small,  since  A,B  are  usually  also  small:  the  coherent 
scatter  signal  component  <&s)>  also  vanishes  for  "geometrical"  reasons. 

3.4  <^>2 

As  can  be  seen  from  (2.25),  we  need  <0^>2  for  the  evaluation  of  detector  performance. 
We  readily  find  from  (3.1 1),  (3.12)  here  that  now 


<0(s)>2=  jB^t’on)2  |  So(t'on)in  I2 

•  { 1  -  cos  2  [CQot'on  —  ^inO'on)  ~  koDo(t'on)] ) 


(3.13) 


where 

B^(t'on)  =  [loko VgOF R0S  (2a02)Fi(t,0n)AvDWRT(2koCXol)e-RG/2]  (3.13a) 

=  {4S0OozkogTgR-4REF^o^  Fl(t'on)AvDWRi(2k0aol)e"RG/2 

.  e"  ((2aox~bxRx)^  +  (2a0y-bxRy)2/B  )/2  j  ^  j  3^) 

with 

A  ref  =  nfc20  VAB  ,  cf.  (A.2-33),  (A.3-7),  [5],  (3.13c) 

a  reference  area,  used  in  the  definition  of  scattering  cross  sections,  vide  Appendix  III  of  [5].  The 
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prninem  facto*  of  (3. 13a,b)  reveal  the  three  components  which  can  degrade  the  magnitude  of  the 
coherent  scatter  signal:  (1)  doppler ", smear' '  -<F,)J  ■  (ii)  smfoce  roughness  (-e’1^)-  (Ui) 
geometry"  (~wrt). 

signal.^'  “  ^  imP°m"t  10  tha'  ^  ®  »  abWe  "*  W'P“d“'  “f  driving  or  input 


3J  The  Scatter  Signal  Covariance  <  flW  dsi 

ft  V/l 


This  quantity  is  essential  for  the  prediction  of  detector  perfomtance  in  purely  incoherent 
reception,  cf.  (2.26)  above,  where  <6<>>>  =  0.  From  (32b),  [5],  we  have,  generally 


(3.14) 


where  now 


*°(X)S  /  *kfT)crtoKaF ;  *  to(f)  s ISin(f)l2 ,  t=  (n-n')At, 

— ©o 


(3.14a) 

is  the  covariance  of  the  normalized  input  signal  envelope  wavefomt,  s(t-e....),  shifted  to  r(=f-f„) 
or  these  narrow  band  signals,  in  the  manner  of  (29b, c),  (5)  and  (3.3)  above.  Note  that  K^O)*  = 
here.  The  quannty  ts  specified  by  (29d)  et  seq.  of  [5),  now  with  the  modified  doppler  temt 

12oto*vWlo\T 

e  -»  exp  ( iko(2oco«yc  +  2ao\yd  -  r0r v0t+  «oR*v0r)t }  (3 j 5) 

for  our  present  model,  cf.  (3.2). 

This  Mg(t)  becomes  generally  Eq.  (51),  (5)  for  the  postulated  composite  wind-wave 
surfaces  in  force  here.  We  have  explicitly 

M$(t)  =  e-^WW  .  R2fT<cri2kiflQ. 'VDT;^ 

*  [  Hg  J i3(Ar)F2(Ar,x)G  e2^0^1-  d(Ar) 

So 

+  ko  /i3(Ar)Mn(o)(Ar,x)F2(Ar,t)GMs(Ar,T)e2^,4rd(Ar)],  (3.16) 

So 

which  is  canonical  in  the  statistics  of  the  surface  elevation,  fc.  Mn(0),  Ms  are  respeedvely  the 
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second-moment  functions  of  the  wave  "tilt-factors"  and  the  surface  soliton  component,  cf.  (49a)- 
(51).  Here  Hq  has  the  form 


HG  * 


+  +  oi  2  \2 

Qv  oz; 

Ctoz/2 


(3.16a) 


"high-frequency"  case,  RG  »  1,  (3.12c);  =  b2  =  (2c<oz)2,  "low-frequency"  case,  Rc  <  1;  (see 
Sec.  HI,  D(2)  of  [5]  for  details).  We  have,  also, 


(A2-23),  [51:  !3(Ar)  =  J  aRT(ri)  aRT(ri  +  Ar)*dn 

—oo 

-  SSS 

WW 

=  (grgR)2  /l  REF  e-*2"4'-*™-"4 .  e-ik°"-bnt. 


J  (3.16b) 


cf.  (A.2-30),  also  (A.2-18)-(A.2-22),  [5],  specifically  for  the  gaussian  beam  patterns  of  (3.12k) 
and  [5].  The  (partial)  contributions  of  the  "doppler  spread"  are  here 

<e“i2koao*fb't>AvD  =  e-(2a0Kic)\2c2/2-(2a0*od)2k02t2/2  .  c~(t0r<JT)2k^2/2-<1oR^R)2k2t2/2  (3  J?) 


where  the  dependence  on  x  alone  stems  from  the  postulated  stationary  of  the  various  random 
processes  involved. 

In  detail,  we  may  distinguish  the  two  important  cases  of  "high-"  and  "low-frequency" 
operation.  For  the  former  we  have  from  (53),  [5]  specifically  for  in  (3.14)  now 


/.  "High-Frequency"  Cases,  Rq»  1: 


M$(t) = (gm?  a  ref 


(tt2  +  Ot2  +  (X2 
1  OX  ^ 

<Xoz/2 


_ 22 _ e-[(2oox-brRx+ia|Tb^c0x)2/2b|oJx+(2a0y-l>rRy)2/2b^C!^y] 

k^^GxOGy 

k2  e^t2/2  Ng,>_inc«x0)Ws(2cx0ikolx)],  (3.18) 


where 
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Rg 


„2 

aGx,y 


s  (boa^/co)2,  Eq.  (A.2-19),  [5]; 

=  [(2ao*Oc)2  +  (2ao*Od)2  +  ( rorOT)2  +  (  ioR’CTR)2]^’  c{-  (3-17) 

=  Eq.  (A.8-29b),  with  (A.8-26a),  [5];  b0  =  cos0ot  +  cos60r  =  2ccoz,V 


=  ^kj,  <y2t,  (A.8-9a,b)  and  (A.8-18),  [5] 


(3.18a) 


NGS-inC(a°)  =  {Eq-  (4g):  Appendix  IV,  A  J5) }  is  the  "tilt  factor," 
while 


Ws(20ojkolT)  =  Eq.  (54)  of  [5]; 


(3.18b) 


this  is  the  wave  number -time  intensity  spectrum  of  the  soliton  surface  ensemble-,  vide  Appendix  IV, 
B-3,4  also. 

For  the  grazing-angle  cases  (<j>0  <  30:  0ot  >  60°),  with  proper  wind  conditions  [U»  >  0(5 
m/s)J  we  get  the  simpler  form  of  (3.14)  explicitly: 


<0(s)0(s)>  =  Re  {  B^DinKoWin  e“i<0oX  }  , 


(<J)0  30°)  (3.19) 


where  now  the  "correlation  function"  B0  is 


lbi  =  a2!?2  S2  G^k4  (gTgR)2A  ref  N  GS-inc^00^2**000  T°T  11  •  Ws(2<Xo±ko  It). 

(3.19a) 


with  a  doppler  spread 


drift,  deep  current,  and 
platform  motion,  (3.18a) 


gravity-wave  phase  modulation 
Oj>M),  vertical  motion  (3.18a) 


+  Rg)  (3.19b) 

AM  (=  4>M)  by  "tilting" 
of  gravity  waves,  (3.18a) . 


The  doppler  shift  (~  vd)  is  given  by  (3.15)  explicitly.  (From  (3.14a)  we  note  again  that  Ko(0)in  = 
1).  The  general  effect  of  B0(t),  which  is  comparable  in  time  scale  to  Ko(T)in,  is  to  "smear,"  as 
well  as  shift,  the  resulting  signal  spectrum,  and  change  its  level. 
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II.  "Low "Frequency"  Cases,  Rq~<  1: 

For  these  cases  we  use  (55a),  [5],  dropping  the  "dc"  term.  The  result  for  (3.14)  is  now 


<0<*>0§>>  =  Re  {  BW(t)ioKo(x)in  c***  } , 


(3.20) 


cf.  (3.19a),  where  now,  with  (49b)  and  (55c)  of  [5] 

B<s>(t),0  SJ  J2-G(i)e2ikoao^Dt-<^^/2  (gTgR)2A  ^ 

•  {4/1  ref  e^x-hTRx)2^  +  (2a0jrbTRy)2/B  +  (b0k0)4WG(2k0a0iyT) 

+  kjNg>_inc(a0)e'^2/2WS(2aoiko/t) } .  (3.20a) 

The  rather  different  structure  of  B0(x)io  vis-H-vis  B0(x)hi  stems  from  the  behavior  of  F2g  with  Rg, 
cf.  (A.8-6,8)  vs.  2),  Sec.  m,  E,  of  [5].  The  first  term  of  {  }  in  (3.20a)  usually  vanishes,  unless 
?0tox-bTRx  =  0,  etc.,  and  we  are  left  with  the  contributions  of  the  large-scale  surface  waves,  which 
usually  dominate,  except  at  very  small  grazing  angles,  since  0OT,oR  ->  ti/2,  or  b0  ->  0,  while 
Ngs-uic  ^8  aGy  s^n4®oT  >  0,  (4»oT  =  Jt/2).  Again,  B0(x)io  spreads  the  spectrum  of  the  received 
scatter  signal.  Note,  however,  that  both  B0(x)hi  and  B0(x)i0  are  independent  of  the  injected  signal 
waveform  in  these  "narrow-band"  cases. 

4.  Elementary  Models  of  Target  Scatter  Signals,  0(T) 

Target  modeling  from  the  viewpoint  of  scattered  radiation,  like  scattering  from  random  ocean 
surfaces  and  bottoms,  is  also  a  venerable  problem.  Here  we  present  a  short  hierarchy  of 
elementary,  i.e.,  approximate  models,  without  claiming  dominant  originality.  The  main  purposes 
of  this  hierarchy  are  (i)  to  provide  specific  forms  of  target  signals,  0W,  needed  in  our  analysis  of 
detection  performance,  and  (ii)  to  indicate  some  further  directions  for  obtaining  realistic  target 
scatter  waveforms,  which  are  ultimately  required  for  signal  processing  applications  generally. 

Our  present  and  preliminary  hierarchy  of  models  may  be  briefly  summarized: 

Model  I :  Narrow-band  signals  impinging  on  a  continuous,  closed  (3-dimensional)  acoustic 

reflecting  surface,  in  the  far-field:  This  case  is  readily  developed  in  a  formal  way, 
directly  from  Eq.  (3.1a)  et  seq.,  above,  by  some  simple,  formal  modifications  of  the 
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directly  from  Eq.  (3.1a)  et  seq.,  above,  by  some  simple,  formal  modifications  of  the 
TSSF,  (3.1b).  See  Sec.  4.1  ff.  and  Fig.  4.1.  Evaluation  of  the  integrals, 
however,  presents  technical  problems. 

Model  la :  This  is  the  same  formulation  as  Model  I  but  with  the  3 -dimensional  surface  replaced 

by  a  suitably  oriented  plane  scattering  surface,  in  the  manner  of  Fig.  4. 1  a. 

Model  II:  Here  the  continuous,  3-dimensional  scattering  surface  is  replaced  by  a  set  of  fixed, 
but  possible  randomly  located,  reradiating  elements  ("emitting  sensors"  in  an  array). 
The  same  narrow-band,  far-field  conditions  are  assumed  here.  This  is  the  "discrete 
facet"  or  "aligned  point  scatterer"  model. 

Model  III:  This  is  Model  I,  but  now  considered  in  the  Fresnel  zone,  as  far  as  phase  variations  in 
the  incident  and  scattered  radiation  is  concerned.  In  this  model  we  can  account  for 
"large"  distributed  targets  (~Lmax)  when  the  Fraunhofer  (far-field)  conditions  Lmax/R 
«  1, 27tLniaxCos6oTAo  «  RoT.oR  (cf.  Eq.  19. 15  of  (23])  are  no  longer  obeyed, 
i.e.,  the  so-called  "bow-aspect." 

Model  Ilia:  This  is  the  partial  near  field  or  Fresnel  case  of  Model  la,  the  plane  surface  version  of 
Model  III. 

Model  IV:  This  is  the  discrete-element  version  of  Model  HI. 

Model  Class  V:  These  involve  broad-band  signals  impinging  on  continuous  or  discrete  reflectors, 
in  the  manner  of  Models  I,  II  above,  in  the  far  field. 

Model  Class  VI:  This  class  is  the  same  as  that  of  Model  I,  but  now  in  the  semi  near  field,  or 

Fresnel  phase  region.  Extensions  include  the  near  field  effects  in  the  amplitude  as 
well. 

All  of  the  models  are  based,  to  varying  degrees  of  approximation,  on  the  basic  boundary  value 
approach  of  classical  theory  [22],  [23],  of  which  the  Kirchoff  or  Tangent  Plane  method  is  here 
used  primarily  [5],  for  frequencies  0(£  0.2  Hz),  cf.  remarks  in  Appendix  IV,  C,  [5]. 

In  addition,  because  the  objects  in  question  are  often  compliant  structures,  which  can 
resonate  when  excited  by  incident  radiation,  an  additional  re-radiation  source,  with  an  associated 
directional  beam  pattern  determined  in  detail  empirically  for  the  type  of  structure  involved,  is  added 
to  the  models  above,  cf.  Sec.  4.4,  Eqs.  (4.23),  (4.24). 

Finally,  when  we  are  concerned  with  the  passive  situation,  only  the  target  as  possible  emitter 
needs  to  be  considered:  this  is  the  well-known  "transient"  mode,  which  can  usually  be  ignored  in 
our  present  concern  with  the  detection  of  active,  reflected  sources  here.  It  should  be  noted, 
nevertheless,  that  there  are  two  main  ways  to  proceed  with  "transient"  signals:  (1)  the  simplest  is  to 
process  the  (partially)  random  class  of  signal  waveforms  as  seen  at  the  receiver,  and  (2)  the 
considerably  more  difficult  one,  namely,  to  proceed  along  the  lines  of  the  present  analysis  (cf.  [5]), 
which  require.,  the  explicit  modeling  of  both  the  source  and  the  various  media  through  and  in 
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which  the  original  signal  is  propagated. 


4.1  Target  Signal— Model  I 

This  model  represents  a  direct,  formal  application  of  our  previous  result  (3.4),  starting  with 
(3.1a).  Now,  instead  of  a  more  or  less  flat,  open  two-dimensional  surface  S0,  cf.  Fig.  3.1,  as 
scattering  element,  we  replace  it  with  a  closed,  three-dimensional  surface  S3,  representing  our  dis¬ 
tributed  target  body.  Equation  (3.4)  is  formally  modified  for  the  narrow-band  signals  used  here  to 


0(?l„.b.  =Re  {SoSoOpn'lin  ^)(iko,t'pn,...)peiC0otPn'  }, 
f.f. 

tpn'  —  nAt  —6  -  Top; 


(4.1) 


where  now,  however,  is  given  by  [cf.  (3.1b)] 

=  iko ^/G^-^^W.tpn)  J(/?  o5)p  (ft.^)  OR-Krp’lfo)  e~ikorp,,2a(oP)  dS3(rp'),  (4.1a) 

£3 

in  which  the  integrand  is  to  be  evaluated  over  the  body's  surface  S3.  Here  we  have  for  the  doppler 
component  (from  Eq.  (24),  [17]) 

D^s  (2ocp*vc)(tpn'+T0p)  +  (2ap»vc)(tpn’+  R0t/c0) 

+  {-OoTP’V0T)tpn'  +  ('ioRp*v0R)(tpn’  +  T0p) }  ,  (4.2a) 

and 

rp'  =  r  +  £p;  r  =  !rr,  Cp  =  »r  *  »p  =  0,  (4.2b) 

so  that  r  and  £p  are  always  perpendicular.  The  directional  wavefront  normals  (2cto->)  2ap  are  now 
2ctop  =  (ioT  -  'oR)p  -  Eq.  (A.  1-3)  of  [5],  where  RqT  -»  RoTp.  RoR  ->  RoRp.  etc.,  by  obvious 
modification  of  the  geometry  of  Fig.  (3.1),  for  the  ocean  wave  surface  vis-i-vis  source  and  receiver, 
ft,  as  before,  is  an  inwardly  directed  normal  to  the  surface  S3.  Similarly,  T0p  =  (RoTp  +  RoRp)/co  is 
the  path  time  from  Or  -»  Op  ->  Or;  also  vp  now  is  the  (constant)  velocity  of  the  body,  again 
referred  to  the  fixed  reference  system  0(x,y,z)T,p,R,  with  (lvl/c0)2  «  1,  [17].  The  geometry  of  the 
body  system  is  shown  in  Fig.  4.1.  The  quantity  Grt,  as  before,  is  the  projection  of  the  beam 
pattern  on  the  surface  S3;  is  the  usual  spreading  factor,  given  by  (A.  1-3),  with  absorption. 
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Fig.  4.1:  Geometry  of  the  scattering  body,  vis-k-vis  the  source  (T)  and  receiver  (R),  cf.  Fig.  3.1. 


The  critical  analytical  problem  here  is  the  evaluation  of  the  surface  integral  (I  S3).  This  is 
usually  very  difficult,  so  that  approximations  or  numerical  methods  must  be  resorted  to.  In 
addition,  the  direction  of  r  (=  y),  namely,  the  orientation  of  the  body,  is  not  usually  known,  and 
must  therefore  be  regarded  as  random  at  the  receiver.  The  quantity  (#o£)p  (ft»2ot{g))  Ort  can  be 
regarded  as  the  effective  beam  pattern  projection  on  the  now  three-dimensional  surface 
analogous  to  the  corresponding  situation  in  the  problem  of  scattering  from  the  wind-wave  surface 
discussed  in  Section  3,  cf.  (3.12j).  For  most  bodies  /?<$  will  be  unity  (as  distinct  from  water  — » 
air  — > water  scattering)  cf.  [23],  pp.  220  et  seq. 

4.1  A:  Coherent  Moments  of 

The  various  first-  and  second-order  moments  of  8$,  needed  for  optimum  threshold  coherent 
detection  algorithms  and  performance,  cf.  Sec.  (2.4),  can  be  written  down  at  once  from  the  above 
(4.1),  (4.2),  in  (3.11),  (3.13),  namely. 


<e(?>  =  Re  {SoSo< W)in  <  F^Wtpn)>peiwo,P«}, 
tjjn’  *  nAt  -e  -  Top; 


(4.3) 


where  the  averages  are  now 
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<  >  =  <  >, 


(4.3a) 


'’aO’®r’  ’ 

cf.  eq.  (3.10).  The  physicogeometric  "modulator,"  <F^>,  of  the  original  signal  becomes 

<f£)j>  =  iko  R  0p  e^V  (Fi)Avd  <WRT(2a0pko)p>ir  (4.4a) 

where/?  0p  =  1,  while  S  p  =  1  (for  the  surface  of  the  body  not  "in  shadow,"  as  seen  by  R,T).  For 
the  doppler  we  use  (4.2a)  in  (3.12d— i).  The  scattering  contribution  of  the  body  itself  is  given  by 
the  effective  wave  number  scatter  spectrum 

<wRT(2appk0)p>ir  =  J <(A-2aop)  QRT(rp'lfo)  e_ik°rP',2.2oP>  .dl3(rp'),  (4.4b) 

where  the  average  is  over  the  body's  orientation.  This  is  the  quantity,  of  course,  which  presents 
the  principal  analytical  difficulties  in  evaluation  for  this  model  and  which  prompts  us  to  seek 
alternative,  but  valid,  approximate  forms,  cf.  Section  4.2  ff. 

The  analogue  of  <0(^>2,  (3.13),  here  is  directly 


(4.5) 


where  now 

Bo(t’on)p  =  'Soko  (^o5)p  Fi(t’on)AvD<  wRr(2appk0)p>,  (4.5a) 

cf.  (3.13a,b).  Again,  "doppler  smear"  (~Fi-avd)  will  degrade  (and  possibly  destroy)  the  coherent 
structure  of  the  body  scatter.  Unlike  the  wave  surface  scatter,  however,  there  is  no  "roughness" 
degradation  here,  usually.  Finally,  note  the  different  path  delays  (~  T0p,  etc.)  and  directionalities 
(~  otop),  etc.:  geometry,  as  always,  plays  an  important  role. 
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4.1-2:  The  Incoherent  Moments  <0^  6$  >;  <6^ $  > 

From  (29),  [5]  we  can  write  the  analogue  of  (3.14)  here  as 


<9W>  =  Si  G^k2  Re  {  Ko(x)i„  Mrt(t)|j  }  . 


(4.6) 


Now  (cf.  (29d)  [5])  and  the  doppler  of  (4.2a),  with  (3.12d)-(3.12i),  (3.17)  suitably  modified 
gives  us 


M$(t)  =  e~i2k°£°*''Dt  •  R?  <(e'2iko2o*rD^>AvD  < WRT(2aopkolO)p>fr  (4.6a) 


in  which 

<WRT(2cXopkolO)p>ir  s  // <(aRT(r’ip  lfo)aRT(r’ip  +  Ar’  lf0)*(  fti-2oc0p)(A2«2a0p)>fr 

£ 

•  e2Lk°^P'Ar'  dXi(r'ip)  dKAr')  (4.6b) 

is  the  wave  number  (intensity)  spectrum  of  the  scattering  body,  analogous  to  (4.4b).  Here  hi  - 
Ai(r’ip),  A2  =  ft2(r’ip  +  Ar"),  and  Ar’  always  connects  two  points  on  the  fixed  surface  X:  Ar’  is 
not  random,  although  r’,  r  are,  in  direction.  Again,  a  central  problem  is  to  evaluate  WRT-p, 

(4.6b).  Thus,  Equation  (4.6)  can  be  expressed  in  the  same  form  as  (3.19)  and  (3.20),  viz., 


<0(t)0(T)>  iRe{  B^T)Ko(T)in  e  ic0°”  }  , 


(4.7) 


where  the  "modulation  effects"  of  the  doppler  shift  and  spread,  and  in  particular  the  scaling  by  the 
scatter  body  (~  Wrt-P),  are  described  by 


B(£\t)  =  a£  G^l  e-2^?*  R?  <c2iko2pP*DT>AvD  < WRT(2aop  I0)p  >tf .  (4.7a) 


Specifically,  we  have  for  the  doppler  spread  here 
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(4.8) 


<e2ikoaop«VDT>AvD  =  e^x2/2 . 

cl  m  [(2aop  •  ac)2  +  (2ao3  *  aPp)2  +  ‘oTp  •  oj)2  +  i0Rp  •  OR)2]k2. 

from  (3. 17),  where  we  have  replaced  gd  by  cpp,  the  vector  spread  for  the  target's  motion. 

For  and  ^0  $8  n^»  we  proceed  in  a  similar  way,  combining  (3.4)  and 

(4.1)  and  averaging  appropriately,  according  to  (4.3a).  The  result  is,  finally, 

<eW>  -  Ctoz^k2  y(GW$R0SR0v  • 

Re{Ko(t  +  Aop)in  e-ic0o(t+AoP)<wRT(2koaol)s  wRT(2koa0p)p>  •  <e_ik°D  oiW#2>AvD}  , 

Aop  =  T0  —  Top,  (4.9) 

with  D(*\  D<§>  given  respectively  by  (3.2),  (4.2a),  and  with  Aop  — >  -Aop  and  — >  -D^, 

->  -D^)  in  (4.9)  for  <0(*)0(*)>.  Here  we  have  as  a  typical  average  for  the  doppler  factor*  in 
(4.9): 


<e-ikoao*V-Avc‘l  +  ^o^P’Ve+^VcOj  +t)> 

c 

=  e“iko^c^lao~(tl+x)a0pl .  Fi(-kotiOo  +  ko(ti+x)0Cop)Avc  (4.9a) 

where  generally  do  *  ctop,  etc.  Thus,  when  ti  — >  tj  +  to,  and  to  — » «»  for  cw  signals,  not  only  is 
(4.9a)  representative  of  a  nonstationary  process  (since  dependent  explicitly  on  ti),  but  Fi  ->  0,  in 
the  usual  cases  where  Avc  *  0.  Again,  as  expected,  doppler  spread  destroys  coherence, 
completely  for  very  long  signals,  as  in  the  case  of  <0(s)>,  <0(T)>,  (3.11),  (4.3)  above,  also.  Even 
for  comparatively  short  signals  doppler  spread  due  to  platform  and  scatter  motion  can  essentially 
destroy  coherence  in  these  far-field  situations,  because  of  the  long  travel  times  involved,  i.e.,  large 
RoT/co,  RoR/co.  T0,  Top,  etc.,  cf.  (3.12f-i). 

For  most  practiced  purposes  [unless  possibly  ctop  =  (Xq  and  T0p  =  T0,  which  occurs  when 
the  target  is  sufficiently  close  to  the  surface  scatter  domain,  cf.  Figs.  1.1  a,b],  the  cross-moments 
<0^  >,  <C  & $  dfi  >  effectively  vanish.  However,  even  when  the  target  (P)  and  surface 

scatter  region  (S0)  are  close,  so  that  <Xop  =  cto,  Top  =  T0  are  good  approximations,  the  doppler  shift 


"'The  epoch  (e)  is  omitted  in  there  slowly  varying  terms,  since  e  is  very  small,  0("RF”  cycle):  the  doppler  shift  in 
the  envelope  S0.i„  is  thus  negligible. 


31 


and  spread  terms  in  (4.9)  become  explicitly 

<e-ik0Doi)+ik0D?2  >Avd  =  e^c  *  “o  +  VrVT4,  'oR‘voR* 

.  <e2iko°o  ‘  f^p  +  4vp)(tl+r+ToR)  -  (  v  d  +  Avd)(t1-ToR)]>Avd^vp  ^  (4 .9b) 

so  that  unless  vp  =  vd  and  Avp  =  Avd,  the  doppler  spread  again  destroys  coherence  between  the 
surface  scatter  signal  and  that  scattered  from  the  target  (p).  Our  conclusion  here,  accordingly,  is 
that/or  most  practical  purposes  the  cross  signal  moments  between  surface  scatter  and  the  target 
vanish,*  e.g., 


<0  W>  =  <e«e$>  *  o. 


(4.10) 


This  very  considerably  simplifies  the  performance  parameter  (CTn)inc  in  (2.26). 

42  Model  la:  An  Approximate  Model  1: 

The  principal  approximation  here  is  the  replacement  of  the  three-dimensional  closed  surface 
I3  of  the  scattering  body  of  Model  I  above  (Sec.  4. 1)  by  an  "equivalent"  flat,  2-dimensional 
scattering  surface,  in  the  manner  of  Fig.  4.2.  Then,  it  is  convenient  to  use  the  coordinates  of  the 
flat  target  surface  to  evaluate  WR-p.p,  with  S  -»S  =  1.  These  coordinates,  in  turn,  need  to  be 
related  to  the  primary  (parallel)  system  of  Ot  and  Or  (cf.  Fig.  4.1  and  Fig.  4.1a).  This  is  done 
with  the  help  of  the  familiar  rotational  transformations  ([24],  Sec.  8): 


h 

tv 
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'x' 
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fy 

mi 

m2 
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03 

*x’  =  y  1 + y  2 + y  3 

fx  =  ix’/ 1  +  •y'ml+  «z*ni 
etc.,  (4.11) 


*This  result  will  not  hold,  however,  in  the  near  field,  with  signals  of  short  duration,  as  then  the  random  doppler 
spread  will  not  have  sufficient  time  to  destroy  the  coherence  between  the  scattering  surface  and  the  target. 
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where  l\,  lj, ...  are  the  direction  cosines  of  Ix\  the  unit  vector,  ix’,  etc.,  in  the  primed  system  of 
the  "flat"  target.  Specifically  (cf.  Fig.  4.  la)  we  have  again  (4.2b) 

r’p  =  r  +  Cp  =  iff  +  ipcp  ;  «r  ■  «p  =  0; 


("inward"  drawn  normal); 

(4.12b) 

and  in  terms  of  the  primary  (x,y,z)  systems,  we  get  direcdy 

f  ‘P  =  ix'  =  »x  cos  <{>p  sin  9p  +  iy  sin  $p  sin  9p  +  iz  cos  =  ixli  +  iyl2  +  Ub,  etc. 

(4.13)  <  ir  =  iy-  =  ix  cos  <j>r  sin  9r  +  iy  sin  <>r  sin  0p  +  lz  cos  0r  , 

l  *r  X  ir  =  iz  =  ix  {sin  <J>p  sin  0p  cos  0r  -  sin  <J)r  sin  0r  cos  0p) 

+  iy  {cos  <j)r  sin  0rcos  0p  -  cos  <J)p  sin  0p  cos  0r) 

+  iz  {cos  <(>p  sin  0p  sin  0r  -  sin  <J>p  sin  0p  cos  0r  sin  0r}  , 

where  (0r,4>r).  (9p.<t>p)  are  the  respective  polar  and  azimuthal  angles  of  r  and  £ p  in  the  primary 
system,  cf.  Fig.  4.1a.  Comparing  the  coefficients  ofix,  fy,  tz  with  the  relations  of  (4.1 1)  shows 
at  once  that  they  are  respectively  the  direction  cosine  //,  etc. 

From  the  above  it  is  clear  that  wrt-P  (4.4b)  can  be  evaluated,  once  #rt  (r'plf0)  is  specified, 
since  n  =  iz,  so  that  (20op)Z'  is  readily  found  to  be 

(2a0p)z'  •  (i0pT  -iopR)  =  iz'  •  L(«x.  iy.  iz)  (4. 14a) 

where  is  i0pT  -ioPR  given  by  (A.  1-3),  [5].  Thus,  we  have  with  the  help  of  (4.1 1),  in  the  bistatic 
cases. 
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y 


(20Cop)z'  =  [{(1  +  RoTp/RopR)  cos4>rp  sin0rp  -  (Lq/RoPr)  cos<|>opL}ni 


+  {(1  +  RoTp/RopR)  sin<j>Tp  sin8jp  -  (l^RopR)  simfopUnz  +  (cos0OTp  +  cos0opR>n3], 

(4.14b) 

with  RopR  given  by  (A.l-2c),  [5],  (with  appropriate  indexing),  where  ni,  n2,  n3  are  the 
coefficients  of  ix,  ly,  tz  in  iz\  (4.13)  above.  The  wave  number  scatter  spectrum  WRj.p,  (4.4b) 
now  becomes 


<WRT((2app)z-]p>jr=  /  (2<Xop)z’  J aRT(r’  lf0)  ^‘V’^a'op-l  dx'dy'  V  { 

Ap  */®o’ V 


(4.15) 


where  2'op-l  =  U'OXoplx’ +  »y'(°k>p)y'.  in  which  the  components  of  (Xop  are  readily  found  from 
(4.1 1),  (4.13),  and  the  general  forms  (A.l-3),  [5],  as  done  above  for  (2cx0p)zi,  (4.14b). 

Next,  we  use  (3.12k)  for  the  specific  gaussian  beam  patterns  of  our  current  models,  with  r 
->  r'  (=  r'p  here)  =  1xx'  +  Iyy',  cf.  (4.12a),  and  (Bjr,  bjR)  -»  (B'tr,  b'jR)  for  the  particular 
geometry  of  T  — » (5  -»  R  now.  All  the  results  for  <0(s)>  in  Section  3  above  cany  over  directly  in 
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form,  with  (2003)2'  replacing  (20o)z',  etc.,  and  with  appropriate  modifications  of  the  doppler 
relations,  geometry,  etc.  See  (3.12/),  (3.13a,b,c),  etc.  Thus,  for  this  simplified  model  la  we  have 
approximate,  explicit  relations  for  the  various  target  moments,  also  with  appropriate  modifications: 


<0W>  =  Eq.  (3.11),  with  (3.12),  [Eqs.  (3.12a,b,c)  omitted],  (3.12j-m);  (4.16a) 

<0(t)>2  =  Eq.  (3.13),  with  (3.13a, b,c);  (4.16b) 

<0$,0^>  =  Eq.  (4.7),  with  <Wrt-3>,  (4.6b),  (4.7a),  replaced  by  <Wrt-3>,  (4.16c) 

(4.17)  below,  which  now  for  gaussian  beam  patterns  is  found 
to  be  (without  the  average  over  ir) 


WRT(2k0oto3)3  =  (2a0p)z'  J ^RT(r'ilfo)  &RT(r'i  +  Ar’lf0)*  e2iko(a’op)±  •  Ar'dr’id(Ar’),  (4.17) 

Ap 

specifically,  cf.  Fig.  (4.1a), 

=  (2<Xop)z'  Il[(Oop)l  -bRT/2,  (Oop)l-bTR/2],  Eq.  (A.2-25),  [5]  (4. 17a) 

=  (2a0p)z'  4(gTgR)2(A'REF)2  ekoVoPU  •  *TR  •  (a'opU  (A.2-16),  [5]  (4.17b) 

with  A'ref  =  =  A'l/2’  and  (a'op)l  =  (aop)i  -b'RT/2.  (4.17c) 


Similarly,  we  have  for  (4.4b)  here 


WRT[2ko«op]p  =  (2otop)z' 


Ap 


ORi(r'lf0)  e2,kor*‘  (a'op)ldr' 


A  A 

=  Eq.  (3.12/, m)  Btr  -»  B'tr,  b'jR,  etc.,  in  the  gaussian  cases.  (4.18) 

With  these  modifications,  then,  we  get  for  the  cross-averages  <0^,0^>,  etc.,  Eq.  (4.9)  with 
appropriate  doppler  adjustments  in  (4.9a,b)  and  Wrt-P  given  by  (4. 1 8): 


<0^),0^>  =  Eq.  (4.9),  with  (4.18),  etc.,  remarks  above.  (4.19) 


Once  again,  via  the  arguments  of  Section  4. 1-2  for  Eqs.  (4.9a),  (4.9b),  we  may  expect  these 
cross-averages  to  vanish,  e.g.,  (4.10)  holds  here,  also,  for  this  simplified  model,  with 
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considerable  simplification  of  the  performance  parameter  on^7n2c  (2.26). 

Finally,  for  Model  la  here  we  also  assume,  in  addition  to  the  far-field  conditions,  that  Ap, 
Fig.  (4.  la),  is  very  large  compared  to  the  wavelength  of  the  incident  radiation,  so  that  "edge 
effects"  may  be  ignored.  (But  see  our  remarks  in  Sec.  4.3  ff.) 


43  Model  II:  "Point-Scatter  Model" 

Here  we  replace  the  continuous  plane  surface  Ap  of  the  scattering  body  by  a  series  of 
scattering  facets  (also  in  the  plane  Ap),  whose  locations  are  variable,  in  the  manner  of  Fig.  4.1b 
below.  This  is  the  "point-scatter"  version  of  Model  la,  whose  practical  valuation  depends  on  the 
empirical  values  we  may  assign  to  the  various  reflection  coefficients  Rom  associated  with  each 
"facet"  about  point  r'mp. 

Again,  we  can  write 


Cm  -  ix'CmJ  rm  -  tyTmi  A  -  iz';  (rm  •  Cm  -  0);  ®  *  2(Xop  -  (20Cop)z’ .  (4.20) 

Here  F^,  Eqs.  (4.1a),  (4.4a),  apply  formally,  with  wRT-p,  (4.4b),  therein  now  replaced  by 

M 

Wj^(20opko)discrete  =  S  (2otop)z'  om  ClRTd-’mp  lf0)  e'^P  *  2£°P,  (4.21) 

m 


where 


Z 


Fig.  4.1b:  Point  Scatter  Target  Model 


r  mp  -  ^x'Cm  +  tylm  (4.21a) 

r  mp  #  2ctop  *  (2otop)x’  Cm  +  (2®op)y'  rm>  (4.21b) 

with  (4.13)  specifying  the  components  iX’,  iy', 
etc.  vis-zi-vis  the  basic  coordinate  system  (x,y,z) 
here,  as  before.  Thus,  is  again  an 

effective  wave  number  scatter  spectrum  of  the 
target  body  P  in  this  model,  which  allows  us  to 
assign  various  weights  /?0m  (and  /?om).  as  well 
as  various  locations  along  the  line  y\  to  each 
"point"-facet,  at  r'mp. 

In  a  similar  way  we  see  that  for  the 
second -order  moments 
<e(x)0(ns)>t  etc.,  cf.  (4.6),  we  have  now 
(5^-4  1): 
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Wj^(2koC(op)discrete  =  <Tz'  *  2ctop)^  X  R  om  R  om'  ClRTd-Vnp  lf0;aRT<r'mp  lfo)*e2‘k°a'0P*(r’m‘-r'm) 

(4.22) 


mm 


for  (4.6a),  (4.6b),  which  may  be  subject  to  the  further  average  <  >ir=sy-  as  noted  above,  when  . 
target  orientation  is  unknown  at  the  receiver. 

The  "point-scatter  model"  as  presented  here  is  in  part  necessarily  phenomenological,  i.e.,  we 
must  assign  reasonable  values  to  such  quantities  as  R om,  Rom  R  om-,  rm,  On.  based  on  empirical 
studies.  Nevertheless,  it  has  the  attractive  features  of  (comparative)  acoustic  simplicity  and 
identification  with  empirical  data.  Also,  "comer  reflector"  or  "high-light"  components  may  be 
included,  as  additional  (random)  scattering  points  in  Ap.  We  add  at  once  that  all  this  is  not  a  new 
concept.  Our  treatment  here,  however,  incorporates  a  number  of  features  not  always  included, 
namely,  doppler  effects  and  bi-static  configurations  with  the  body  p  plane,  along  with  source  and 
receiver  beam  patterns. 


4.4.  A  Resonance  Component 

In  addition  to  the  scattered  radiation  there  may  often  be  a  "resonance"  component,  whereby 
the  scattering  body  reemits  radiation  in  a  number  of  "normal"  or  resonant  modes,  when  the  incident 
radiation  is  sufficiently  intense  and  consists  of  short  pulses.  Here  we  model  the  response  simply 
as  a  response  function,  Y0(s)p,  cf.  Sec.  2.2-5, 920,  [7],  associated  with  an  indigenous  beam 
pattern  Qp(vp  lfQ).  Then,  this  added  component,  appearing  in  (4.1a)  in  Models  I,  la,  II,  is 
simply 

Y0(s)pap(v0p)  {/?0pV G(p)aRT(r'p)e'soTo'(so/co>^PT  "  1opR)‘r'P‘(so/co)D^)5(r'p-0)  }  , 

(4.23) 

with  (vGp)  =  rx'fo/c0,  cf.  (4.2a)  for  the  doppler,  where  Y0p  is  the  narrow-band  resonance  response, 
at  f0,  of  the  target,  if  any,  and  where  Op  is  determined  by  the  target  structure.  This  latter  quantity 
must  generally  be  determined  empirically  and  is,  of  course,  a  function  of  orientation  (fx’,  (4.13)), 
which  is  usually  random  vis-si-vis  the  observer. 

Accordingly,  we  must  add  to  WRT-p,  and  WRT-p,  (4.4b),  (4.6b)  et  seq.,  the  quantity  (4.23). 
We  have,  then,  for  Model  I: 


Eq.  (4.4b),  <  tfT)  >:  <wrt-P>{j.  -»  {  <wRT-p>  +  Yo(s)pap(vop)QRT(0  lf0)  } 


Eq.  (4.6b),  <  ( <WRT.p>^-» 


(4.24a) 
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|<Wrt_(j  >+2Re{<wRT_jj  >Y0_pQp  +YopClpflRx(0)2  }|  (4.2b) 

Eq.  (4.9),  <&snl&y>:  <WRT.p>lr -»  Eq.  (4.24a)  (4.24c) 

Similarly,  WRx.jj  for  Models  la,  II,  are  modified,  cf.  (4.15),  (4.18);  (4.17);  (4.21),  (4.22),  by 
adding  the  appropriate  "resonance"  component,  cf.  (4.24a-c). 

5.  The  "Telecommunication  Model":  T  ->  x;  T  -»  S-4t: 

In  what  we  call  here  the  "telecommunication"  model,  because  of  its  analogy  to  the 
communication  situations  in  which  signals  are  transmitted  directly  to  the  receiver  -  here  on  the 
scattering  body  (x)  -  transmission  proceeds  as  indicated:  (i)  directly  T  -»  t,  or  (ii)  off  the  ocean 
surface  (S),  to  the  receiver,  e.g.,  T  -» S  x,  in  the  manner  of  Fig.  1.1b  above. 

The  detection  formulation  of  Secs.  2.3, 2.4  apply  here,  provided: 

A Case  I,T->  x:  0«  ->  6;  0(*>  ->  0,  in  (2.19)-(2.26); 

(5.1)  I  Case  II,  T  -*S  —>  x:  0W  — >  0;  0(s)  *  0,  with  receiver  R  replaced  by  the 
\  body  (P),  e.g.  Or  -» Op,  RoR  -»  Ro^r,  etc.;  in  the 

scatter  results  of  Section  3,  cf.  Fig.  3.1,  with 
corresponding  modifications  of  the  doppler,  (3.2)  et  seq. 

Hj:  0  +  0^s)  +  Namb  vs.  Ho:  Namb; 

0W  ->8+  0(s) in  (2.13a);  add  <0^,0(f»  in  (2.13b); 

add  <0(s)>2  in  (2. 19); 

v  add  <0^,0^>2  in  (2.1 1) 

Case  I  is  the  more  common  situation  in  practice,  although  Case  II  can  be  important  if  the  receiver  is 
comparatively  close  to  the  surface,  so  that  the  direct  and  scatter  paths  are  approximately  the  same. 
Then  the  receiver  can  use  the  scatter,  or  multipath  component  to  enhance  detection: 

6.  Results,  Remarks,  and  Next  Steps 

Here  we  briefly  summarize  the  principal  results,  with  comments,  and  indicate  some 
possible  next  steps  in  the  analysis: 

6.1  Principal  Results 

In  the  preceding  we  have  determined 

I.  The  parameters  of  the  optimum  threshold  signal  detectors  [Secs.  2.1, 2.2],  and 

II.  The  associated  parameters  for  optimum  threshold  performance  [Secs.  2.3, 2.4],  where 
general  nongaussian  (plus  gaussian)  noise  is  present  and  when  signal-dependent  noise  is  present, 
produced  here  primarily  by  scattering  off  the  ocean  surface  [Sec.  3].  The  results  are  canonical, 
i.e.,  independent  of  the  specific  statistics  of  the  ambient  (i.e.,  background)  noise,  and 
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are  similarly  canonical  in  signal  structure. 

These  results  are  also  quite  general:  "bistatic"  as  well  as  "monostatic"  geometries  are 
included.  Particular  attention  is  given  to  the  role  of  doppler  in  the  platforms  and  scatterers. 

III.  Target  modeling  is  explicitly  treated  (Sec.  4):  rather  general  models  are  outlined,  and 
approximated  by  effective  two-dimensional  (flat)  scattering  surfaces  and  point-scatter 
elements  in  a  plane.  Ihese  models  are  quasi-phenomenological,  requiring  in  practice 
calibration  with  empirical  data. 

Throughout  we  have  emphasized  realistic  structures  and  have  based  much  of  the  detailed 

results  on  the  recent  analysis  presented  in  [5],  which  should  be  used  in  conjunction  with  the 

present  work  for  the  fullest  utility.  It  should  be  noted,  moreover,  that  in  all  the  cases  treated  here 

* 

the  input  signal  (envelope)  S0(t)in.  is  explicitly  a  factor  in  the  parameters  describing  the  threshold 
algorithms  and  performance,  either  directly,  or  squared,  or  in  the  input  signal  covariance,  cf.  (3.4), 
(3.8a),  (3.11),  (3.13),  (3.14a),  (3.19),  (3.20),  for  the  signal-dependent  scatter,  and  for  the  target 
returns,  cf.  (4.1),  (4.3),  (4.5),  (4.6),  (4.7),  (4.9).  This  explicit  and  not  too  surprising  separation 
of  the  input  signal  from  the  physical  and  geometric  factors  which  embody  the  scattering 
mechanisms  and  their  geometries,  doppler,  etc.,  is  strictly  the  consequence  of  the  narrow  band 
assumption,  in  conjunction  with  the  far-field  condition,  assumed  throughout  this  initial  effort. 
When  cither  or  both  of  these  conditions  are  removed,  i.e.,  broad-band  signals  and/or  Fresnel 
geometries  are  invoked,  then  the  present  analysis  must  be  appropriately  extended:  the  input  signal 
no  longer  is  factorable  from  the  physico-geometric  factors,  but  appears  as  a  linear  functional,  either 
directly  as  such  or  as  a  functional  of  the  signal  covariance. 

6.2  Next  Steps 

Various  next  steps  need  to  be  considered.  Among  them  are  (not  necessarily  in  the  order 
indicated): 

A.  Applications  to  LPI problems,  with  particular  attention  to  the  role  of  physico-geometric 
factors  involved  in  the  wave  surface  and  target  scattering  situations; 

B .  Extension  of  the  analysis  to  the  Fresnel  region,  in  order  to  account  for  targets  in  the 
beam  orientation,  where  the  Fraunhofer  condition  breaks  down; 

C.  Calculation  of  representative  results,  for  typical  geometries; 

D.  Quantitative  study  of  the  doppler  effects-, 

E.  Formal  extension  of  the  analysis  to  broad  band  signals,  as  a  framework  for  numerical 
analysis; 

F.  Application  of  the  current  results  to  the  estimation  of  waveform  and  target  level. 

Other  important  extensions  and  applications  will  present  themselves  in  the  course  of  the 
development  of  A-F  above. 
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Appendix  A.l  Definitions  and  Descriptions 

The  principal  purpose  of  this  Appendix  is  to  provide  definitions  of  the  various  terms  and 
elements  of  Eqs.  (3.1a,b),  and  others,  which  appear  in  the  text  (Many  of  these  may  also  be  found 
in  Appendix  I  of  [5],  as  well  as  in  Parts  I— IH  of  [6].) 


We  begin  with  Eqs.  (3.1a,b): 


(A.l-1) 

(A.  1-2) 

e 

(A.  1-3) 

GO) 

(A.l-4) 

T0 

(A.l-5) 

Co 

(A.  1-6) 

Wo 

(A.  1-7) 

Ro 

(A. 1-8) 

S 

(A.  1-9) 

So 

(A.l-10) 

A  -  ftG 

(A.l-1 1) 

foT>  «oR 

(A.l-12) 

20q 

(A.l-13) 

VD 

(A. 1-14) 

f(  )fr 
J  2m 

Brl 

(A.l-15) 

Vt,  Vr 

(A.  1-16) 

{ — )g+s 

(A. 1-17) 

C(r,t) 

(A. 1-18) 

?oPT,  «oPR 

=  TSSF  =The  'Total  Surface  Spreading  Function,"  Eq.  (3.1b),  which 
describes  the  effects  of  the  random  moving  wave  surface  £(r,t)  on  the 
incident  and  scattered  acoustic  signals 
=  the  signal  "epoch,"  as  measured  at  the  receiver 
=  e‘2ao^oT°/(47t)4R0jR0‘R,  the  "geometric"  or  spreading  factor, 
including  the  effects  of  absorption  for  wave  surface  scatter 
=  (Rot  +  RqR)/c0  =  path  delay  Ot  -»  S(O')  -»  Or,  cf.  Fig.  1,  [5] 

=  (av.)  speed  of  propagation  in  the  water  medium 
=  2rtf0  =  carrier  signal  angular  frequency 
=  plane  wave  reflection  coefficient  for  the  air-water  interface 
=  shadowing  function,  cf.  Appendix  IV,  [5] 

=  A  =  portion  of  the  mean  ocean  surface  jointly  "illuminated"  by  the 
transmitting  and  receiving  beams,  Ot,  Or 
=  unit  normal  to  the  gravity-capillary  wave  surface,  e.g.,  A  =  (ix£x  + 
tyCy  -  tz)nz,  C  =  CGI  Cf.  Eq.  (23),  [5],  with  nz  =  (1  +  g  +  ^-1/2 
=  unit  vectors  along  R0t,  R0R,  cf.  Fig.  (la,b)  here  and  Fig.  1,  [5] 

=  foT  -  *oR.  cf.  Appendix  I,  [5] 

=  a  doppler  velocity,  which  can  include  platform  motion,  surface  drift, 

etc.;  cf.  Sec.  VIII,  Appendixes  VII,  Vm  of  [5] 

i»+d(>0) 

=  a  Bromwich  contour  (of  integration)  =  f  (  )  — ,  which 

-i~+d  2m 

allows  transient  or  finite  duration  signal  inputs,  as  well  as  those  of 
infinite  duration  (d  -» 0) 

=  volumes  occupied  by  the  transmitting  and  receiving  arrays,  cf.  Fig. 
3.1. 

=  quantity  [in  (3.1b),  (3.6b)],  associated  with  the  composite  surface  G 
(=  gravity  +  capillary  waves)  and  S  (=  postulated  surface  soliton 
ensemble  [5],  [6] 

=  wave  surface  elevations 

=  unit  vectors  along  R<>pT,  RoflR.  cf.  Figs.  4.1,  4.1a 
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(A.l-19)  Fg,  Fg 

(A.  1-20)  G(p5 
(A.l-21) 

(A.  1-22)  RG,g 


=  TSSFs  for  target  body  (£3),  and  two-dimensional  (flat)  target  (£2), 
cf.  (A.l-1). 

_  e-2a(1^CoToP/(47t)4Ro|rRo^R,  "geometric  factor"  for  target,  cf.  (A.  1-3) 
=  doppler  terms,  cf.  (3.2),  (4.2a),  for  wave  surface  and  target 
=  Rayleigh  numbers 
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